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Analyses  are  made  of  the  mutual  interactions  between  shock  structure 
and  the  sidewall  laminar  boundary- layer  and  their  effects  on  the  quasi- steady 
flat-plate  laminar  boundary  layer  in  ionizing  argon  shock-tube  flows.  The 
mutual  interactions  are  studied  using  effective  quasi- one-dimensional  equations 
derived  from  an  area-averaged- flow  concept  in  a finite-area  shock  tube.  The 
effects  of  mass,  momentum  and  energy  non-uniformities  and  the  wall  dissipations 
in  the  ionization  and  relaxation  regions  on  the  argon  shock  structure  are  dis- 
cussed. The  new  results  obtained  for  shock  structure,  shock-tube  laminar  side- 
wall  and  quasi- steady  flat-plate  boundary-layer  flows  are  coopered  with  dual- 
wavelength interferometric  data  obtained  from  the  OTIAS  10  cm  x 18  cm  Hypervelocity 
Shock  Tube.  It  is  shown  that  the  difference  between  the  results  obtained  from  the 
present  method  and  those  obtained  by  Enomoto  based  on  Mirels'  perfect-gas 
boundary-layer  solutions  are  significant  for  lower  shock  Mach  nunbers  (Ms  ~ 13) 
where  the  relaxation  lengths  are  large  (~  10  cm) . In  general,  the  present 
results  agree  better  with  our  experimental  data  than  our  previous  results  for 
uncoupled  ionizing  flows. 
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Notation 


shock-tube  cross-sectional  area 

(Vi w 
P6P6 

local  skin  friction  coefficient 

specific  heat  at  constant  pressure 
kA 

-j-  hydraulic  diameter 

anbipolar-dif fusion  coefficient 

plasma  total  enthalpy 

Boltzmann  constant 

shock-tube  perimeter 

shock-wave  Mach  nunfcer 

atomic  mass  of  argon 

local  Nussalt  number 

electron-nuirber  density 

electron-nunber  density  production-rate 

pressure 

elastic  energy-transfer  rate  to  electrons 
inelastic  energy- transfer  rate  to  electrons 
radiation  energy  Iobs 
shock-tube  sidewall  heat  conduction  rate 
plasma  heat  conduction  energy  flux 
plasma  diffusive-energy  flux 
electron  heat- conduction  energy- flux 
electron  diffusive-energy  flux 


v 


R gas  constant 

Re  Reynolds  number 

T temperature 

Tj  characteristic  ionization  temperature 

t time 

u velocity  component  in  x- direction 

V particle  velocity 

v velocity  component  in  y-direction 

x distance  along  shock-tube  sidewall 

y distance  normal  to  shock-tube  sidewall 


a 


h 

P 

T 

Subscripts 

a 

e 

i 


w 

o 

6 

Operator 

<f> 


degree  of  ionization 
specific  heat  ratio 

boundary  layer  thicknesses;  see  Eqs.  (28a)~(28j) 

ion-diffusive  velocity 

thermal-conduction  coefficient 

mixture  viscosity  coefficient 

total-plasma  density 

shear  stress 


atom 

electron 

ion 

conditions  at  the  shock-tube  sidewall 
conditions  in  front  of  shock  wave 
condition  at  the  boundary  layer  edge 


dA  area-averaged  values 
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1 . INTRODUCTION 


In  our  previous  analyses  of  ionizing  shock-wave  structures  (Glass 
et  al.  Refs.  1,  2)  and  the  ionizing  boundary-layer  flows  induced  by  a strong 
shock  wave  (Liu  et  al,  Refs.  3>  4) , the  coupled  effects  between  inviscid  and 
viscous  flows  in  a shock  tube  were  neglected.  As  shown  by  Glass  and  Patterson 
(Ref.  5)  and  Mirels  (Ref.  6),  the  flow  between  the  shock  wave  and  the  contact 
surface  in  an  actual  shock  tube  is  nonuniform  owing  to  ..he  growth  of  the  shock- 
tube  sic'ewall  boundary  layer.  The  gross  features  of  shock-wave  structure  in  an 
ionizing  gas  are  affected  by  the  sidewall  boundary  layer.  Consequently,  the 
induced  ionizing  boundary  layers  on  a flat  plate  in  the  so-called  quasi-uniform 
flows  axe  affected  by  these  mutual  interactions.  This  nonuniformity  has  to  be 
considered  when  interpreting  shock- tube  data  in  aerodynamic  or  chemical-kinetic 
studies. 


The  in^i or t ant  effects  of  the  growth  of  the  sidewall  boundary  layer 
on  the  freestream  .flow  are  the  induced  wall  shearing  stress,  heat  transfers  and 
the  consequent  nonuniformities  in  the  flow.  In  order  to  take  into  account  the 
flow  nonuniformities  in  a shock-tube,  Mirels  (Refs.  6,  7,  8)  has  obtained  corre- 
lation formulae  for  a perfect  gas,  "hich  are  often  used  by  many  researchers  in 
gasdynamics.  These  relations  are  based  on  similarity  solutions  of  the  boundary- 
layer  equations  treating  the  boundary  layer  as  an  aerodynamic  sink.  Local  - 
similarity  assumptions  cannot  be  applied  to  those  cases  where  the  variations  of 
the  freestream- flow  quantities  are  important  and  cannot  be  neglected.  As  shown 
by  Liu  et  al  (Refs.  3,  4)  the  similarity  assumption  is  not  valid  for  the 
electron-tenqperature  and  electron-mmber-density  profiles  for  the  sidewall- 
boundary-layer flow  where  ionizing  nonequilibrium  phenomena  occur  in  the 
freestream  flow  behind  the  shock  front  and  the  variations  of  the  freestream 
conditions  for  the  sidewall  boundary  layer  are  significant.  Also,  only  mass 
conservation  was  considered  by  Mirels  in  the  study  of  flow  nonuniformity. 

The  effects  of  the  sidewall  boundary  layer  on  the  freestream  momentum  and  energy 
equations  were  neglected.  The  validity  of  Mirels'  correlation  formulae  cease 
to  apply  to  an  ionizing-gas  flow,  although  they  were  used  by  Enomoto  (Ref . 9) > 
McLaren  and  Hobson  (Ref.  10)  and  Brabbs  and  Belles  (Ref.  11)  in  their  studies 
of  shock  structure  in  an  ionizing  gas.  Enomoto  studied  the  sidewall  boundary- 
layer  effects  on  ionizing  shock  structure  in  argon.  He  applied  Mirels'  (Ref.  6) 
perfect-gas-correlation  formulae  for  flow  nonuniformity  of  the  mass  mixture. 

Flow  nonuniformities  in  the  momentum  and  energy  equations  were  entirely  neglected 
without  giving  any  reasons.  Nevertheless,  he  found  the  important  physical  result 
that  the  relaxation  length  is  significantly  reduced  in  tubes  of  decreasing  cross- 
sectional  area  owing  to  the  sidewall-boundary- layer  growth. 

The  sidewall  boundary- layer  growth  generates  both  compression  and 
rarefaction  wavelets.  However  the  net  effect  is  that  of  a rarefaction  wave 
which  attenuates  the  shock  front;  see  Trimpl  and  Cohen  (Ref.  12).  Flow  uni- 
formity in  ionizing  gases  in  the  relaxation  regions  induced  by  incident  and 
reflected  shock  waves  is  still  far  from  being  completely  understood..  In 
principle,  the  analytical  approach  can  be  made  by  taking  into  account  the  inter- 
action between  the  boundary-layer  growth  and  the  development  of  a two  or  three- 
dimensional  inviscid- flow  model  as  well  as  oy  taking  into  account  their  inter- 
action with  unsteady  flow  effects  such  as  shock-wave  attenuation  or  contact 
surface  acceleration. 

For  example,  Hubbard  and  de  Boer  (Ref.  13)  studied  two-dimensional 
flow  nonuniformities  in  a perfect  gas  based  on  the  basic  assumption  that  the 
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boundary- layer  parameters  do  not  change  very  much  in  the  region  where  two-dimen- 
sional effects  are  of  importance.  They  showed  that  the  flow  in  the  inviscid 
region  is  nearly  one-dimensional  at  a distance  from  the  shock  front  greater 
than  the  shock-tube  radius.  Demmig  (Ref.  l4)  studied  one-dimensional  ionization 
relaxation  in  krypton  under  influences  of  weakly  attenuating  shock  fronts  and 
sidewall  boundary-layer  effects.  He  calculated  the  flow  field  behind  an  experi- 
mentally determined  attenuating  shock  front  by  taking  into  account  only  the 
boundary-layer  displacement  thickness  growth  in  a way  similar  to  Enomoto's  method. 
He  showed  the  combined  effects  of  boundary- layer  growth  and  shock-wave  attenuation 
on  the  ionizing  krypton  shock-wave  structure  at  a shock  Mach  number  of  about  10. 

In  the  usual  one-dimensional  model,  the  boundary  layer  is  assumed  thin 
compared  to  the  shock-tube  radius  and  is  replaced  by  a sink  distribution.  The 
inviscid  core  is  assumed  to  be  quasi-one-dimensional  in  the  sense  that  flow 
variations  are  assumed  to  occur  only  in  the  streamwise  direction.  However, 
as  Glass  et  al  (Ref.  2)  showed  experimentally  that  the  electron- cascade  front 
approaches  the  translational  shock  front  near  the  sidewall,  the  mutual  inter- 
action between  ionization  mechanisms  and  the  flow  field  is  not  simple. 

In  the  present  analysis,  the  inviscid  core  field  is  assumed  to  be 
effectively  one-dimensional  by  using  an  area- averaging  process  instead  of 
treating  the  boundary  layer  as  a sink  distribution  as  for  the  perfect-gas  case. 

The  purpose  of  this  work  is  to  extend  our  previous  studies  of  shock- 
wave  structure  and  boundary-layer  flows  for  ionizing  argon  by  including  the 
mutual  interactions  between  the  inviscid  freestream  and  viscous  sidewall 
boundary-layer  flows.  The  present  study  of  the  mutual  interactions  is  based 
on  the  effective  one-dimensional-flow  equation  derived  from  the  flow  area- 
averaged  concept,  where  the  flow  nonuniformities  in  the  freestream  mass, 
momentum  and  energy  equations  and  the  wall  dissipation  terms  are  considered. 

A finite-difference  scheme  is  applied  to  the  ionizing  boundary- layer  equations. 

An  examination  is  made  of  the  role  of  the  sidewall-boundary-layer  growth  on  the 
behaviour  of  the  shock- structure  ionization  and  radiation- cooling  regions. 

The  results  of  the  quasi- steady  flat-plate  and  the  sidewall  boundary- 
layer  flow  induced  by  a strong  shock  wave  in  a finite  area  shock  tube  are 
re-examined  and  compared  with  our  dual-wavelength  interferometric  boundary- layer 
data  obtained  from  the  UEIAS  10  cm  x 18  cm  Hypervelocity  Shock  Tube.  It  is  shown 
that  the  difference  between  the  results  obtained  from  the  present  method  and  the 
simplified  Mirels'  one-dimensional  equations  are  even  more  significant  at  lower 
shock  Mach  number  where  the  relaxation  lengths  are  larger.  In  general,  the 
present  analysis  gives  much  better  agreement  with  our  experimental  data  by 
taking  into  account  the  mutual  interactions  between  the  inviscid  and  viscous 
flows. 


2 . THEORETICAL  CONSIDERATIONS 


2.1  Basic  Equations 

Two  methods  have  been  widely  applied  in  obtaining  the  effective 
quasi-one-dimensional- flow  equations  in  gas  mixtures:  (l)  control  volume 
method,  and  (2)  flow  area-averaged  method.  In  the  control  volume  method,  the 
integral  form  of  the  conservation  law  is  applied  to  the  mixture  control  volume. 
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The  effective  equations  are  obtained  by  introducing  definitions  of  area- averaged 
properties  and  limiting  the  length  of  the  control  volume  to  zero.  In  the  flow 
averaged  method,  local  instantaneous  conservation  equations  are  obtained  from 
the  law  of  conservation.  The  effective  quasi-one-dimensional-flow  equations 
are  obtained  by  applying  an  area- average  on  the  local  instantaneous  conservation 
equations . Identical  results  are  obtained  from  both  methods  for  a homogeneous 
gas-mixture  flow. 

In  this  work,  the  effective  quasi -one- dimensional- flow  equations 
for  an  ionized  gas  in  a shock  tube  are  directly  obtained  by  applying  the  flow 
area-average  on  the  local  instantaneous  macroscopic  balance  equations.  Consider 
a gas  mixture  made  up  of  atcms,  singly- ionized  ions  and  electrons.  The  macro- 
scopic balance  equations  for  the  plasma  mixture  are  given  by  Appleton  and  Bray 
(Ref.  16)  as  follows: 


Mass  conservation: 


d 

3t 


P + v • (Pf)  = 


(1) 


Plasma  momentum  equation: 


pV*  + V • (pv^  + VP  = V • ((iW) 


(2) 


Plasma  energy  conservation: 

pH  + V • (pvfe)  = - V • (qc  + q^)  + V(V*  • ^ (3) 

Electron  mass  conservation: 

^ ap  + V • (apV^  = m&ne  + V • (4) 

Electron  energy  conservation: 

apCT^  + V • (pV&My  - (V*  * 5V^)7P€ 

- 7 ' <■%.  * + «.!  * <W 

where  the  definitions  of  the  variables  and  other  terms  are  the  same  as  those  in 
Refs.  1 and  2. 

Consider  a mixture  in  a constant  area  of  a shock  tube  with  cross- 
sectional  area  A.  Integration  of  Eqs.  (l)-(5)  over  the  flow  area  A at  any 


f 


instant  in  time  gives: 


J ||  dA  + J V-  (pV*)cLA  = 0 


(6) 


/It  pV'dA  + J V-(pV^  + p)dA  =J  V-(mW*)dA  (7) 

A A A 

J ^ pHdA  + J V(pHV*)dA  = - J V(qc  + q^)dA  + J V(V%W)dA  - J Q^dA 

A A A A A (8) 

J'  ^ apdA  + J'  V-  (apV*)dA  = maJ^  ngdA  + J'  V*(pDmWt)dA  (9) 

A A A A 

J (apCpTe)dA+  J V*  (pV&CpTe  )dA  = J (V*  + SV^-^dA 

A 

-/ V'(?oe  + 4)"  + / <«el  + <10) 


where  the  conduction-  and  diffusive- energy  flux  are  given  by: 


10  - - + 


and 


^d 

q, 


ce 


(RT!  + CpTe)p^ 


AW 
e e 


(u) 


*de  ' - V^m® 

Using  the  special  forms  of  the  Leibnitz  and  Gauss  theorems  and 
talcing  the  dot  product  to  get  the  x- component  yields: 


Mass  conservation: 


<p>  + <pu>  = ° 


(12) 
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Plasma  momentum  equation: 


<pu  + p>  = - <T>) 

hr'° 


(13) 


Plasma  energy  conservation: 


It  «*>  + & <wao  - o>,  - «"»„  - -V 


(14) 


Electron  mass  conservation: 


It  <P0£>  + |E  <apu>  = ma  <"e>  - <pDm  ^ >w 


(15) 


Electron  energy  conservation: 

It  <ocCPTe>  + 5 <pu*pV  - <(u  + 6ui>  sr  > 

+ <«el  + W * <*.>» 


(16) 


where  the  operator  <f>  = i J fdA  is  introduced,  and  the  boundary- layer  approxi- 

A 

mation  with  neglecting  terms  of  0(52)  is  made. 

The  wall  dissipation  terms  are  defined  as  follows: 


<T>*  ■ s ( “ 5 X 
V» 

<pDm  Ij  >»  = I ( If  X 


(17) 


<UT> 


^e'w-  A ^ce  + qde^w 
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Equations  (12) -(l6)  are  the  basic  equations  for  the  unsteady  area-averaged  flow 
model.  Demmig  (Ref.  1^)  studied  the  effects  of  boundary  layer  growth  and 
unsteady  shock-wave  motion  on  ionizing  krypton  shock-wave  structure  by  solving 
Eqs.  (12) -(l6).  He  took  into  account  only  the  boundary  layer  displacement 
thickness  growth  in  Eq.  (12),  neglecting  other  boundary  layer  effects  and  the 
wall  dissipation  terms  in  Eqs.  (13) - (l6)  as  well  as  a radiation-energy  loss 
term  in  Eq.  (l4) . He  calculated  the  flow  structure  behind  an  experimentally 
obtained  attenuating  shock  wave  by  an  inverse  method,  so  that  inportant  unsteady 
effects  propagating  across  the  contact  front  were  not  taken  into  account.  In 
his  experiments,  the  shock  wave  Mach  number  varied  from  10.6  to  10.1  over  a 
distance  of  2.27  m in  a 5.2  cm  x 5.2  cm  shock  tube.  On  the  other  hand,  in  the 
UTIAS  10  cm  x 18  cm  Hypervelocity  Shock  Tube,  the  shock  attenuation  rates 
AMs/L  are  about  0.2  (Mach  nuatoer/m)  for  shock  Mach  numbers  Mg  = 13  and  16. 
Therefore,  throughout  the  present  analysis,  it  is  acceptable  that  the  shock 
front  propagates  with  constant  velocity  and  shock-fixed  coordinates  axe 
applicable . 


In  the  shock- fixed  coordinate,  the  quas i - one - dimen  si onal  flow 
equations  are  obtained  from  Eqs.  (l2)-(l6): 


Mass  conservation: 


d_ 

dx 


<f$uf>  = 0 


(18) 


Plasma  momentum  equation: 

h <(^  + p>  “ - J ( * ly  \ (19) 

where  the  pressure,  p,  is  constant  across  the  cross  sectional  area  of  the  shock 
tube. 

Plasma  energy  conservation: 


Degree  of  ionization: 


h «*«>-  <4>„-<v  -± 

S “ “a  - I ( ®m  t?  ) 


Electron  energy  conservation: 


dP, 


BE  <p*Ve>  = <(u  + 'V  3ST  > 

+ <Qel  + %nel>  ‘ <^>w 


(20) 


(21) 


(22) 
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The  sidewall  boundary  layer  effects  are  thus  taken  into  consideration 
in  the  area-averaged  terms  of  the  above  equations.  Equations  (l8)-(22)  finally 
become  (for  details  of  the  derivation  see  Appendix  A): 

(23) 

(23a) 

Plasma  momentum  equation: 

|;pu2(A-  ♦ A Sjg-  -l(n  ^ <*> 

Plasma  energy  conservation: 

I It  pu(A  - 5uL)  + rai  S «“< A - y)  *55  cu3<a  - 6eL> 

■ ^ - VA  - ¥■>  - L (“u^)«  (25) 

Electron  mass  conservation: 

It  P>“CA  - ‘ maie'A  ' 6iL>  - 1 ( f®.  ® )„  (26> 

Electron  energy  conservation: 


Mass  conservation: 

pu(A  - B*L)  « 0 
or 

pu(A  - 6*L)  = p u A 


| R ^ ptt*Ce(A  - Qj,  L)  + p2RTe(A  - 6n  L)  u(A  - 6uL) 
e e 


-«*el"W<A-BSL>-L^.  <27> 

The  various  boundary  layer  thicknesses  defined  and  used  in  Eqs. 

(22) -(27),  which  are  related  to  the  known  definitions  of  the  boundary 
layer  thicknesses  (for  example , Schlichting,  Ref.  18),  are  given  as  follows: 

Displacement  thickness: 

mf  ( 1 ’ ) **  ^28a^ 


where  the  prised.  variables  are  those  across  the  boundary  layer  [see  Eq.  (A3)]: 
Momentum  thickness: 

<28b) 

Velocity  thickness: 

6na‘f(1-f')dy  (28c) 

Degree  of  Ionization  thickness: 

(2M) 

Energy- dissipation  thickness: 

5» '/ C1  ~e^)ay  ‘ “*  ' 7")  ^ <26,> 

Radiation  energy  loss  thickness: 


(28f) 


Reaction  rate  thickness: 


Electron  number  density  thickness: 


v/C1*?)* 


Electron  energy- transfer  rate  thickness: 


Electron  terqperature  thickness: 

p'u'a'T' 


(28g) 


(28h) 


(281) 


r / P'u'a'T'  \ r . / a'T!  \ 

\'jr(1”5afri)dy's*t/^‘  (l-«7)*  (283) 
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As  shown  in  the  shock-tube  sidewall  boundary- layer  analysis  in  ionizing 
argon  by  Liu  et  al  (Refs.  3,  4),  if  the  electron  temperature  is  assumed  constant 
across  the  boundary  layer,  the  electron  temperature  thickness  6re  is  approximately 
equal  to  the  degree  of  ionization  thickness  6a,  that  is, 


6t  ^ 6 
xe  a 

It  should  be  noted  that  in  the  shock- fixed  coordinate  6*,  6^,  &u  and  6e 
are  negative  and  6r»  &n  , &ne  &Q  are  positive,  whereas  the  sign  of  6a  depends 

on  the  type  of  boundary-layer  flow,  that  is,  frozen  or  nonequilibrium.  The  growth 
and  details  of  these  boundary- layer  thicknesses  are  shown  in  Appendix  B.  The 
dissipation  terms  - shearing  stress  at  the  shock-tube  sidewall;  wall  heat  trans- 
fers due  to  the  mixture  of  atoms,  ions  and  electrons,  and  the  energy  deposition  at 
the  wall  due  to  the  diffusion  of  ions  - are  explained  in  Appendix  B. 

It  is  worth  noting  that  the  existence  of  the  electron- cascade  front  in 
the  ionizing  inviscid  core  gives  rise  to  the  non  uni form  boundary-layer  develop- 
ment with  distance  quite  differently  from  the  boundary- layer  development  in 
perfect  argon.  Some  of  the  defined  boundary-layer  thicknesses  are  maxima  or 
minima  at  the  electron- cascade  front  and  eventually  influence  the  shock-wave 
structure . 


Enamoto  (Ref.  9)  first  showed  that  the  relaxation  distance  was  reduced 
due  to  the  existence  of  the  shock- wave  induced  sidewall  boundary  layer.  His 
analytical  model  accounted  for  the  effects  due  only  to  the  displacement  thickness 
of  the  sidewall  boundary  layer.  His  boundary-layer  displacement  thickness  was 
based  on  Mirels’  similar  boundary-layer  analysis  for  perfect  argon  (Ref.  6): 


A - 6*L  = A 1 - 


where  N = 0.5  for  the  laminar  boundary-layer  flow  and  J^n  is  the  separating  distance 
between  an  incident  shock  wave  and  a contact  front  calculated  from  the  formulas 
given  by  Mirels  (Ref.  6).  However,  Enamoto  disregarded  the  other  mutual  inter- 
action terms  such  as  the  other  boundary-layer  thicknesses  and  the  wall  dissipation 
terms  which  appeared  in  the  source  terms  of  the  quasi- one- dimensional  equations, 
Eqs.  (24) -(27).  He  also  neglected  the  radiation  energy  loss  term  Qr  in  the  energy 
conservation  equation. 

The  same  discussions  sure  applicable  to  Demmig's  model  (Ref.  l4) . He 
studied  ionizing  krypton  shock-wave  structures  at  shock  Mach  nuufcer  of  about  10 
and  initial  pressure  of  8 torr.  He  considered  the  effects  of  the  growth  of  a 
sidewall  boundary  layer  as  well  as  effects  of  unsteady  shock-wave  behaviour,  and 
solved  the  unsteady  one-dimensional  flow  equations  by  applying  the  characteristic 
method  to  the  so-called  inverse  problem  with  an  experimentally  given  shock-wave 
trajectory.  However,  even  in  his  boundary-layer  model,  only  the  growth  of  the 
boundary- layer  displacement  thickness  obtained  from  the  similar  boundary- layer 
analysis  for  perfect  krypton  was  considered.  Other  boundary- layer  effects  as 
well  as  wall  dissipation  terms  used  in  the  present  study  are  neglected.  Further, 
although  shock  attenuation  effects  were  considered  in  Demmig's  analysis,  the 
radiation  energy-loss  term  QR,was  neglected  in  the  energy  conservation  equation. 


(O? 


(29) 


9 


The  boundary- layer  displacement  thickness  that  Demmig  obtained  from 
an  analysis  given  by  Dem'yanov  (Ref.  15)  is  written  as 


(29a) 


where  Jusdt  - x is  the  distance  from  the  unsteady  shock  front,  d and  Re  are  the 


diameter  of  the  shock  tube  and  the  Reynolds  number  referring  to  the  invlsbid 
uniform  flow  region  and  d,  respectively.  Dem'yanov  (Ref.  15)  assumed  a velocity 
profile  across  the  boundary  layer  as 


where  ti  -J'  ~ dy  and  corresponds  to  t|  at  the  outer  edge 

He  used  the  above  mentioned  velocity  profile  and  solved  the 
behind  a shock  wave. 


of  the  boundary  layer, 
nonuniformity  flow 


Brabbs  et  al  (Ref.  11)  used  the  same  formula  as  Eq.  (29)  to  demonstrate 
the  effects  of  flow  nonuniformities  induced  by  a laminar  sidewall  boundary  layer 
on  the  chemical  kinetic  study  in  a shock  tube.  Belles  et  al  (Ref.  17)  also 
vised  the  same  formula  for  a turbulent  boundary  layer  by  putting  N = 0.8  in  Eq. 
(29). 


Consequently,  only  the  mass  conservation  equation,  Eq.  (23),  is  conmon 
to  Enomoto's  basic  equations.  In  the  present  study  the  displacement  thickness, 
other  boundary  layer  thicknesses  and  wall  dissipation  terms  were  obtained  by 
solving  the  two-dimensional  ionizing  boundary-layer  equations.  The  profiles 
along  the  shock- tube  sidewall  distance  were  significantly  different  from  those 
resulting  from  a similar  boundary-layer  analysis  such  as  given  by  Mirels  (Ref.  6). 
These  profiles  used  in  this  study  are  shown  in  Appendix  B. 

Details  of  an  implicit  six-point  finite  difference  method  for  solving 
the  two- temperature  nonequilibrium  boundary- layer  equations  in  ionizing  argon 
were  given  by  Liu  et  al  (Refs.  3»  4).  Comparison  between  the  first-order 
results  where  the  mutual  interaction  between  inviscid  shock-wave  structures  and 
boundary- layer  flows  were  not  taken  into  account,  and  the  experimental  data  for 
both  quasi- steady  flat-plate  and  sidewall  boundary- layer  flows  were  also  given. 

Once  all  the  boundary  layer  thicknesses  and  wall  dissipation  tmnns  are 
known  from  the  analysis  of  boundary-layer  flows  in  ionizing  argon,  Eqs. 
(23)-(27)  can  be  numerically  solved.  Glass  et  al  (Refs.  1,  2)  solved  Eqs. 
(23)-(27)  for  shock-wave  structures  in  ionizing  argon  and  krypton  where  contri- 
butions from  the  sidewall  boundatry  layers  and  the  wall  dissipation  terms  were 
not  taken  into  consideration.  The  results  were  compared  with  interferometric 
data  in  some  detail  in  Refs.  1 and  2. 

By  rearranging  Eqs.  (23) -(27)  we  have 
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32? 


The  momentum  conservation  equation  is  given  by 

A - 6_L 


1 dp 

■ns  d$ 

pu 


DL  / 

“( 


d . A - 5»L 
dx  Z A - 6mL 


1 du  \ L ( „ du  'N 


(32) 


The  electron  energy  equation  is  rewritten  as 

1 dTe  2 A - 5UL 
Te  dT  = I A - ' 


A - 5UL  A - &neL 


d_  . A lduV2A~  VQel4  Snel 
dx  A - 6mL  " u dx  ) 3 A - 6^L  puoKTe 


^^={!(oD»^)w*^;(xh^r)w}  (33> 


Since  the  electron  energy-transfer  rate  due  to  inelastic  collisions 
is  larger  than  that  due  to  elastic  collisions  in  the  electron  tenperature  range 
now  under  discussion,  the  following  approximation  can  be  made  across  the  boundary 
layer: 

ftel.  * S-nel  _ ^inel  = ^neU 
^el  + ^inel  ^inel  (n) 

“ C 


since  < Q^nel»  8X1,1  <5el  < ^inel* 
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Consequently,  from  the  definitions  of  the  boundary-layer  thicknesses,  Eqs.  (28g) 
and  (28i) , it  is  possible  to  assume  the  following  relationship: 


6n 


(34) 


Variations  of  both  5q  and  5^  along  the  shock-tube  sidewall  are  compared  for 
Ms  = 13*1  and  15.9  in  Appendix  B,  Fig.  b8.  Finally,  by  introducing  the  above 
approximation,  Eq.  (34),  we  can  rewrite  the  electron  energy- transfer  rate  due 
to  the  inelastic  collision  as  follows  (Glass  et  al,  Ref.  2): 

W-k(”e>6(II+lTe)  <35) 

Equations  (29)-(32)  together  with  the  mass  conservation  equation,  Eq. 
(23a)  and  the  equation  of  state 


p = pR(Ta+OTe)  (36) 

can  be  numerically  solved  with  given  initial  conditions  by  using,  for  instance, 
the  Runge-Kutta-Gill  method.  It  is  worth  noting  that  if  the  mutual  interactions 
between  invlscid  and  viscous  flows  are  not  taken  into  consideration,  the  above- 
mentioned  equations  are  exactly  consistent  with  the  basic  equations  of  Glass  et 
al  (Ref.  2)  . 

2.2  Initial  Conditions 


If  all  the  boundary  layer  properties  and  the  wall  dissipation  terms 
known  from  a shock- tube  sidewall  boundary- layer  analysis,  the  basic  equations 
shock-wave  structures,  Eqs.  (30)-(33)>  (23a)  and  (36),  can  be  solved  for  the 
unknowns  p,  Ta>  Te>  u,  a and  p with  the  following  initial  conditions: 


and 


<7  - 1)M  2 + 2 

S 

( 7 + DM/ 


E_ 

Po 


27Mg2  - ( 7 - 1) 

— 


a * o 

(7  + l)Ms2 
(r  - 1)MS2  + 2 


(37) 


aro 

of 
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where  y = 5/3  and  the  subscript  o denotes  the  condition  in  front  of  the  shock 
wave  in  the  shock- fixed  coordinate. 

2.3  Numerical  Procedures 

The  procedures  of  the  present  iteration  scheme  axe  described  as  follows: 

(1)  Calculate  the  first-order  ionizing  shock-wave  structure  by  neglecting  the 
mutual  interaction  terms  such  as  the  boundary-layer  thicknesses  and  the 
wall  dissipation  terms. 

(2)  Using  these  shock-wave  structure  results,  obtain  all  the  boundary-layer 
thicknesses  and  wall  dissipation  terms  by  solving  the  laminar  boundary- 
layer  equations  in  ionizing  argon. 

(3)  Applying  these  boundary- layer  thicknesses  and  wall  dissipation  terms, 
solve  the  area-averaged  quasi-one-dimensional  shock-wave- structure 
equations . 

(4)  Repeat  procedures  (2)  and  (3)  until  the  iterative  solutions  converge  to 
a given  criterion. 

There  is  no  guarantee  that  the  present  iterative  scheme  will  converge.  However, 
from  our  experience  the  solutions  converge  quite  quickly,  within  two  or  three 
iterations. 

As  an  example,  variations  of  the  boundary-layer  displacement  thickness 
with  distance  x at  each  iteration  nuntoer  are  shown  in  Figs.  1(a)  and  l(b)  for 
shock-wave  Mach  numbers  13.I  and  15*9*  It  is  found  that  the  solutions  almost 
converge  at  the  iteration  numbers  of  2 for  Ms  = 15.9  and  3 for  Mg  = 13.1. 

2.4  Discussions 

In  the  present  study,  an  iterative  scheme  for  solving  Eqs.  (23) -(27) 
together  with  the  boundary-layer  equations  was  proposed.  The  shock-wave  struc- 
tures, sidewall  boundary-layer  flows  and  the  quasi-steady  flat-plate  boundary- 
layer  flows  were  re-examined  and  compared  with  the  experiments  given  in  Refs. 

1~4. 


In  the  summary  of  the  theoretical  considerations,  the  following 
differences  are  pointed  out  between  the  present  study  and  the  previous  works: 

(1)  In  the  previous  works  by  Enomoto  (Ref.  9)  for  ionizing  argon,  Demmig 
(Ref.  l4)  for  ionizing  krypton  and  Brabbs  (Ref.  11)  for  chemical  kinetic  studies, 
the  boundary-layer  effects  are  represented  only  by  the  boundary-layer  displace- 
ment thickness.  The  mutual  interactions  between  the  inviscid  shock-wave  structure 
and  the  boundary-layer  flows,  such  as  various  boundary- layer  thicknesses  as  well 
as  the  wall  dissipation  terms  in  the  momentum  and  energy  conservation  equations 
and  so  on,  aire  considered  in  this  work. 

(2)  Similar  solutions  obtained  by  Mirels  (Ref.  6)  were  applied  by  Enomoto 
and  others  in  the  evaluation  of  the  boundary- layer  displacement  thickness.  In 
the  present  study  the  exact  boundary-layer  equations  were  solved  by  an  implicit 
six-point  finite  difference  scheme.  In  addition,  the  variations  of  transport 
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properties  across  the  boundary  layer  were  taken  into  consideration. 

There  were  no  concurrent  analyses  or  experiments  to  show  that  Mirels' 
formulation  (Ref.  6)  of  the  boundary- layer  effects  was  a proper  basis  for  making 
corrections  to  the  previous  analytical  results  or  experimental  data  for  ionizing 
gases.  Enomoto  made  use  of  Mirels'  boundary-layer  formulation  to  study  the 
ionizing  argon  shock-wave  structure.  There  were  still  large  differences  between 
his  analysis  and  the  experiments. 

3.  EFFECTS  OF  SIDEWALL  BOUNDARY  LAYER  ON  SHOCK-WAVE  STRUCTURE 

3.1  General  Considerations 

In  our  previous  studies  of  ionizing  argon  shock-wave  structure  (Glass 
et  al,  Ref.  2) , it  was  shown  that  the  experiments  for  higher  shock-wave  Mach 
number  (Mg  « l6)  are  in  good  agreement  with  the  analysis.  However  the  agreement 
between  analysis  and  experiments  was  not  satisfactory  for  the  lower  shock  Mach 
nuntoer  case  (Mg  « 13)  comparing  ionization  relaxation  lengths  and  the  electron 
number  density  at  the  electron  cascade  front.  Enamoto  (Ref.  9)  showed  that  due 
to  the  boundary- layer  displacement  thickness  effect,  the  ionizing  relaxation 
length  was  reduced  and,  especially  in  the  lower  shock-wave  Mach  number , reduction 
of  the  relaxation  length  became  significant.  However,  his  postulated  analytical 
model  seems  to  be  physically  oversimplified.  In  this  chapter  the  results  are 
compared  with  the  experimental  data  and  the  previous  analyses. 

3.2  Comparison  of  Numerical  Results  with  Experimental  Results 

Figure  2 shows  comparisons  between  the  various  analyses  and  the  inter- 
ferometric experimental  data  given  by  Glass  et  al  (Refs.  1,  2),  for  the  electron 
number  density  profile  of  the  lower  shock-wave  Mach  numiber,  Mg  13*  The  result 
of  the  first  order  solution  neglecting  the  sidewall  boundary-layer  effects,  is 
shown  by  curve  A.  Curve  B shows  the  result  obtained  by  taking  into  account  only 
the  boundary-layer  displacement  thickness  effects,  calculated  from  the  similar 
boundary- layer  analysis  for  frozen  argon;  this  is  similar  to  Enamoto' s model. 
Coefficients  of  the  growth  of  the  boundary- layer  displacement  thickness  -&*/** 
are  shown  in  Fig.  B2  of  Appendix  B.  The  boundary- layer  displacement  thickness 
used  in  curve  B is  written  as  follows: 


-&*  = 0.1006  ^ for  M = 13.1,  pQ  = 5.16  torr  and  Tq  = 300  K, 

and 

-6*  = 0.1064  1 for  M = 15.1,  p =5.10  torr  and  T = 298  K. 

s o o 

Curve  C is  the  first  iterative  result  of  the  present  method  where  the  boundary 
layer  equations  were  solved  exactly,  and  not  only  the  displacement  thickness  but 
also  other  boundary  layer  effects  and  the  wall  dissipation  terms  were  applied  in 
the  analysis  of  shock-wave  structures.  Curve  D shows  the  higher  order  iterative 
result  of  the  present  method.  It  is  obvious  that  the  difference  between  the 
present  results  (curve  D)  and  Enomoto's  model  (curve  B)  is  significant  for  this 
lower  shock  Mach  number. 

The  relaxation  length,  xjj,  is  drastically  reduced  in  Enomoto's  model, 
xe  = 7.8  cm,  and  the  corresponding  result  for  the  first-order  model  is  about 
10  cm,  whereas  the  present  method  gives  xe  = 8.93  cm.  The  experimental  value 
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is  about  xj;  = 9*0  cm.  The  agreement  between  the  present  analysis  and  the  experi- 
mental data  is  significantly  improved  at  lower  shock-wave  Mach  number  Ms  = 13.1. 

These  results  are  tabulated  in  Table  l(  a)  and  l(b)  and  compared  with  other 
results  with  different  parametric  values. 

The  electron  number  density  at  the  electron- cascade  front  predicted  by 
using  Enomoto's  model  gives  a much  higher  value  than  that  predicted  by  the  present 
analysis.  The  reduction  of  the  relaxation  length  and  the  ove i shoot  of  the  electron 
number  density  in  Enomoto's  model  can  be  explained  as  follows.  In  the  first-order 
analysis  neglecting  the  mutual  interaction  effects  the  mass  flux  pu  is  always 
smaller  than  pQus  as  shown  in  Eq . ( 23a) : 

\ 

pu  < pu  *-L-  = pQUs  (38) 


where  the  boundary- layer  displacement  thickness  (6*  < 0)  acts  as  though  it 
enlarges  the  shock- tube  cross-sectional  area  in  the  shock  fixed  coordinate. 

The  production  of  the  electron  number  density,  described  in  Eq.  (26) , is 
as  follows: 


dn 

m n m n 

e 

a e a e 

dx 

pu  poUs 

(39) 


I Therefore  the  quasi -equilibrium  peak  value  of  the  electron  number  density  is 

roughly  proportional  to  the  inverse  of  the  mass  flux  pu.  In  Enomoto's  model, 
the  relaxation  length  is  always  shortened  and  the  degree  of  ionisation  is  over- 
reached. However,  in  the  present  model  the  overshoot  of  the  degree  of  ionisation 
is  comparatively  suppressed  and  consequently  the  relaxation  length  is  not  so 
drastically  shortened  as  Enomoto's  model  predicts. 

The  comparison  between  the  analyses  and  the  experimental  data  of  the 
total  plasma  density  is  plotted  in  Fig.  3*  The  total  plasma  density  initially 
increases  gradually  from  its  frozen  value,  suddenly  increases  near  the  electron- 
cascade  front,  arid  then  reaches  the  quasi-equilibrium  value  as  shown  by  arrows 
in  Fig.  3-  The  total  plasma  density  gradient  is  a maximum  near  the  electron- 
cascade  front  where  the  electron  number  density  production  rate  ne  is  a maximum. 
Due  to  the  radiation  energy  loss,  the  total  plasma  density  slowly  increases 
behind  the  electron- cascade  front.  This  will  be  explained  later  together  with 
the  interpretation  of  the  pressure  profiles.  It  is  shown  that  much  better 
agreement  between  the  present  model  and  experiments  is  obtained. 


A comparison  of  the  degree  of  ionisation  is  also  shown  in  Fig.  4.  The 
relaxation  length  is  reduced  by  about  12$  due  to  the  sidewall  boundary  layer 
effects  at  this  lower  shock-wave  Mach  number  Ms  = 13. 1.  Arrows  on  the  curves 
and  the  experimental  data  indicate  termination  of  the  relaxation  length.  Good 
agreement  between  the  present  model  and  the  experiment  is  obtained  with  respect 
to  the  relaxation  length.  However,  the  predicted  quasi-equilibrium  values 
determined  from  the  present  model,  curve  D,  are  larger  than  the  experimental 
results  by  about  25$.  The  general  tendency  of  the  data  points  is  similar  to  the 
prediction  and  they  are  almost  parallel  to  the  curve  D at  the  quasi-equilibrium 
region.  The  reason  for  this  discrepancy  between  the  analysis  and  experiments  in 
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the  lower  shock  Mach  number  case  is  not  clear.  However,  as  found  later  for  the 
case  of  higher  shock  Mach  number,  agreement  of  a/ag  between  the  analysis  and 
experiments  is  good.  In  the  present  analysis,  the  most  reasonable  formulae  and 
the  most  reliable  coefficients  to  them  are  used  for  the  excitational  cross-section 
for  atom-atom  and  electron-atom  collisions,  the  momentum  transfer  cross-  sections 
for  electron- atom  and  electron-ion  collisions  and  the  radiation  energy  loss  term 
(Ref.  2).  There  may  be  analytically  undetermined  factors  in  estimating  these 
cross-sections  and  the  radiation  energy  loss  term.  Therefore,  it  is  not  clear  how 
these  unknown  factors  depend  on  the  shock  Mach  number. 

Normalized  pressure  profiles  p/p0>  and  temperature  profiles  are  shown  in 
Figs.  5 and  6,  respectively.  Curves  A,  B and  D correspond  to  the  first-order 
results,  Enomoto's  model,  and  the  present  model,  respectively. 


In  the  first-order  results  the  pressure  initially  increases  gradually 
from  its  frozen  value,  and  then  suddenly  increases  to  its  quasi-equilibrium  value. 
The  pressure  gradient  is  a maximum  near  the  electron- cascade  front  where  the 
electron  number  density  production  rate  ne  is  a maximum.  Over  the  quasi-equili- 
brium region,  the  pressure  still  increases  slowly.  The  reasons  are  explained  in 
the  following  argument. 

If  we  rewrite  Eq.  (31)  by  taking  into  account  only  the  boundary- layer  dis- 
placement thickness  for  the  sake  of  simplicity, 


1 du 
u dx 


RT-, 


d — T mane  Sr 

2 i„( A - 6*L)  ♦ -£■  -JT  ♦ -2j 

u pu 


1.5 


) 


(*♦0) 


where  ^ = pu  /yp  < 1 gives  the  local  flow  Mach  number  behind  the  shock  front 
which  is  less  than  unity  in  the  shock- fixed  coordinate.  The  pressure  gradient 
is  also  written  as 
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As  shown  in  Eq.  (4o) , the  velocity  gradient  is  composed  of  contributions 
from  (l)  the  sidewall  boundary- layer  displacement  thickness  growth  d/dx  in(A  - 5*L) , 
(2)  the  electron  number  density  production  rate,  ne,  and  (3)  the  radiation- energy 
losses  Qr.  Initially  the  contribution  from  the  boundary- layer  displacement  thick- 
ness effects  are  greater.  Near  the  electron-cascade  front  the  contribution  from 
the  electron  number  density  proauction  rate  is  more  prominent.  Finally,  behind 
the  electron- cascade  front  the  contribution  from  the  radiation- energy  loss  is 
significant . 


Since  the  flow  field  behind  the  shock  front  in  the  shock-fixed  coordinates 
is  subsonic,  and  the  boundary- layer  displacement  thickness  acts  to  expand  the  cross- 
sectional  area  of  a shock  tube,  the  sign  of  velocity  gradient*  Eq.  (4o)>  is  always 
negative.  Physically  this  indicates  that  due  to  subsonic  flow  behind  the  shock 
wave  in  the  shock- fixed  coordinate,  the  above  mentioned  effects  encourage  the 


velocity  profile  Us  - u to  keep  on  decreasing  with  distance  x.  In  other  words, 
the  particle  velocity  u is  increasing  with  distance  from  the  shock  front  in  the 
laboratory- fixed  coordinate  (see  Glass  et  al,  Ref.  2).  Therefore,  as  obviously 
seen  from  the  momentum  equation,  Eq.  (4l) , the  pressure  should  keep  on  increasing 
with  distance  behind  the  shock  front.  The  curves  B and  D in  Fig.  5 show  compara- 
tively sharp  initial  pressure  rises  from  the  frozen  value  as  a result  of  taking 
into  account  the  boundary-layer  growth  effects.  The  present  results  (curve  D) 
lie  between  the  first-order  results  (curve  A)  and  Enomoto's  model  (curve  B) 
except  at  the  beginning. 

The  same  interpretation  is  applicable  to  the  total  plasma  density  profile. 
The  total  plasma  density  is  inversely  proportional  to  the  velocity;  see  Eq.  (23a) . 
Consequently,  the  atom  temperature  Tn  and  the  electron  temperature  Te  are  over- 
reached in  curve  D as  shown  in  Fig.  6,  and  are  very  much  over-reached  in  Enomoto's 
model,  curve  B.  Atom  and  electron  temperature-overshoots  subsequently  give  rise 
to  more  atom-atom  collisions  and  as  a result  the  electron  number  density  in  front 
of  the  electron  cascade  front  increases.  Therefore  the  relaxation  length  is 
significantly  reduced.  It  is  found  that  Enomoto's  model  gives  an  overestimation 
to  the  reduction  of  the  relaxation  length.  This  interpretation  is  of  course 
consistent  with  the  previous  one  that  the  reduction  of  the  relaxation  length  can 
be  explained  from  the  point  of  view  of  boundary  layer  displacement  thickness 
growth  and  the  reduction  of  the  mass  flux. 

It  is  worth  noting  that  an  initial  sharp  pressure  rise  and  a temperature 
overshoot  are  peculiar  to  the  analysis  where  the  growth  of  boundary- layer  dis- 
placement thickness  is  taken  into  account.  Demmig  (Ref.  l4)  studied  combined 
effects  of  the  sidewall  boundary-layer  displacement  thickness  and  the  unsteady 
shock-wave  motion  on  ionizing  krypton  shock  wave  structure.  He  showed  that  the 
temperature  overshoot  occurs  despite  shock-wave  attenuation. 

One  of  the  interesting  features  of  ionizing  argon  shock  wave  structures 
is  the  existence  of  an  incubation  period  in  which  a gas  particle  passes  from  the 
translational  shock  front  to  the  electron- cascade  front.  In  the  incubation  period, 
no  significant  changes  in  the  flow  properties  occur  but  internal  electronic  states 
are  excited  and  the  electron  number  density  is  very  slightly  increasing.  After 
this  period  all  the  flow  properties  change  very  suddenly.  The  electron- cascade 
front  appears.  Therefore  it  is  physically  acceptable  that  any  change  in  the 

dynamic  flow  properties  such  as  a pressure  rise  induced  by  the  boundary-layer 

growth  would  sharply  influence  the  whole  shock  wave  structure  and  eventually 
reduce  the  relaxation  length. 

Experimental  data  obtained  from  shock  tubes  with  various  configurations 
are  compared  by  means  of  using  the  hydraulic  diameter  D.  Here  we  compared  our 
predicted  values  with  experimental  data  obtained  by  the  IFTIAS  10  cm  x l8  cm 
Hypervelocity  Shock  Tube,  so  that  we  employ  the  standard  D = 4a/L  = 12.86  cm 
for  a 10  cm  x 18  cm  shock  tube. 

In  order  to  examine  the  shock  tube  diameter  effects  on  the  ionizing  shock 
wave  structure,  computations  were  done  for  these  different  diameter  shock  tubes: 

Do. 5 = 6.43  cm,  D1.0  = 12.86  cm,  Dp .0  = 25.73  cm  and  D ->  00  . The  computational 

results  for  the  infinite  shock  tube  diameter  Da,  is  equivalent  to  those  of  the 
first  order  analysis  since,  as  shown  in  Eqs.  (23)-(27),  the  terms  representing 
the  mutual  interaction  effects  are  inversely  proportional  to  the  hydraulic 
diameter.  On  the  other  hand,  in  a small  diameter  shock  tube  the  mutual  inter- 
action effects  are  even  greater,  and  in  the  extreme  case  of  a very  small  diameter 
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shock  tube,  the  area-averaged  quasi -one-dimensional  analysis  is  no  longer  valid. 

In  this  case,  a two  or  three-dimensional  flow  model  and  the  subsequent  strong 
shock  attenuation  effect  should  be  considered.  Sometimes  a very  small  diameter 
shock  tube  is  used  as  a common  shock  tube  or  a detonation  tube.  For  example, 
Enomoto  (Ref.  9)  compared  the  reduction  of  relaxation  length  from  a one-inch 
diameter  shock  tube . However,  it  is  not  clear  whether  the  present  method 
can  be  applied  to  this  small  diameter  shock  tube  or  not.  In  Table  1(a)  the 

quasi-equilibrium  peak  value  of  Qfg,  ngg,  pg  and  the  relaxation  length  xg  are 

tabulated  at  Mg  = 13 .1,  p0  = 5.l6  torr  and  T0  = 300  K for  various  shock-tube 
diameters.  In  addition  both  nonequilibrium  and  frozen  boundary-layer  profiles 

are  compared.  However  it  is  known  from  the  results  of  the  shock-tube  sidewall 

boundary  layer  that  the  nonequilibrium  profile  is  closer  to  the  experimental 
data  ^sse  Chapter  4 of  this  report  and  Liu  and  Glass,  Ref.  4).  Therefore,  the 
higher-order  iteration  is  applied  only  in  the  nonequilibrium  flow. 

Figure  7 shows  the  result  of  the  electron  number  density  for  the  various 
hydraulic  diameters  where  rigg  is  the  equilibrium  electron  number  density  for 
= 12.86  cm.  As  shown  in  Fig.  7,  the  relaxation  length  decreases  and  the 
electron  number  density  at  the  electron- cascade  front  increases  as  the  shock- 
tube  diameter  is  made  smaller. 

As  a second  example  at  higher  shock- wave  Mach  number  Ms  = 15 .9 > the 
experimental  and  analytical  results  for  the  electron  number  density  are  shown 
in  Fig.  8.  In  this  case  the  difference  between  the  present  model  (as  shown 
in  curve  D)  and  Enomoto' s model  (curve  B)  is  smaller  than  that  in  the  previous 
low  shock-wave  Mach  number  case  Ms  = 13 .1 • The  reason  is  taat  effects  of  the 
various  boundary  layer  thicknesses  and  wall  dissipation  terms  on  the  shock-wave 
structure  are  not  so  important  since  the  quasi-equilibrium  temperature  is  higher 
and  the  relaxation  length  is  shorter,  xg  = 2 cm  for  the  higher  shock  Mach  number 

Mg  = 15-9. 


In  Fig.  9 the  total-plasma  density  profiles  are  shown.  Good  agreement 
is  obtained  between  the  present  model  and  the  experiments.  Results  for  the 
degree  of  ionization  profile  are  shown  in  Fig.  10  for  conparison.  Pressure  and 
temperature  profiles  are  shown  in  Figs.  11  and  12,  respectively.  The  same 
remarks  as  for  the  lower  shock  Mach  number  are  applied  to  this  case.  The  quasi- 
equilibrium peak  values  axe  shown  in  Table  l(b)  for  various  shock-tube  diameters. 
The  effects  of  the  shock-tube  diameter  on  the  shock-wave  structure  are  plotted 
in  Fig.  13.  As  expected,  effects  of  the  shock-tube  diameter  on  the  shock-wave 
structure  are  smaller,  since  the  relaxation  length  is  shorter  in  this  case  so 
that  the  t:'  me  duration  in  which  the  boundary  layer  influences  the  structure  is 
shorter.  It  is  found  that  the  profiles  for  a D2.0  shock  tube  and  an  infinite 
diameter  shock  tube  are  almost  the  same.  The  profile  for  a = 12.86  cm  shock 

tube  is  very  close  to  that  for  Doo.  This  shows  that  the  17TIAS  HM>ervelocity 
Shock  Tube  gives  almost  one-dimensional  shock-wave  structures  at  higher  shock 
Mach  numbers.  In  the  first-order  analysis,  the  relaxation  length  is  2.15  cm, 
whereas  in  the  present  model  xg  = 2.03  cm  after  the  iteration  number  of  2.  Conse- 
quently, at  this  higher  shock  wave  Mach  number,  the  boundary- layer  effect  reduces 
the  relaxation  length  by  only  6%  whereas  it  is  about  12%  at  the  previous  lower 
shock  wave  Mach  number. 

3.3  Discussions 


Comparing  our  area-averaged  quasi -one-dimensional  model  with  Enomoto' s 
model  and  the  first-order  analysis  as  well  as  interferometric  data,  the  following 
discussions  are  pointed  out: 
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(1)  The  formerly  selected  value  of  argon-argon  collisional  cross-section 
constant  of  10“19  eV/cm£  (Glass  et  al,  Ref.  2),  which  was  used  throughout  the 
present  study,  is  accurate.  The  argon-argon  collisional  cross-section  constant 
is  known  to  affect  the  relaxation  length.  Goal  agreement  was  obtained  between 
experiments  and  the  present  analysis  for  both  Ms  = 13. 1 and  15-9. 

(2)  The  relaxation  length  is  significantly  reduced  due  to  the  sidewall 
boundary-layer  growth  as  predicted  by  Enomoto  (Ref.  9)-  The  reduction  of  the 
relaxation  length  is  more  pronounced  for  the  lower  shock-wave  Mach  number 

(Ms  = 13*1)  where  the  length  is  much  longer  (xe  w 9.0  cm)  than  that  (xe  = 2.0  cm) 
for  the  higher  shock-wave  Mach  number  (Ms  = 15  *9)  • 

(3)  The  model  suggested  by  Enomoto  over-predicted  the  effect  of  boundary- 
layer  growth  on  the  shock-wave  structure . Therefore  the  application  of  Mirels* 
method  to  ionizing  gas  flows  is  inadequate,  except  for  a rough  estimation  of  the 
effects  of  the  sidewall  boundary  layer  on  the  shock-wave  structure. 

(4)  The  effect  of  the  sidewall  boundary  layer  on  the  shock-wave  structure  is 
more  pronounced  in  a shock  tube  of  smaller  diameter.  For  higher  shock-wave 
Mach  numbers  the  reduction  of  the  relaxation  length  is  small  enough  for  us  to 
conclude  that  the  UTIAS  10  cm  x 18  cm  Hypervelocity  Shock  Tube, with  which  the 
present  experiments  were  conducted,  is  large  enough  for  us  to  ignore  the  effect 
of  the  boundary  layer  on  experiments  in  ionizing  argon  flows. 


4.  SHOCK-TUBE  SIDEWALL  BOUNDARY -LAYER  FLOWS 
4 .1  General  Considerations 


Shock-tube  sidewall  boundary-layer  flows  in  ionizing  argon,  without 
taking  into  account  the  mutual  interactions  between  the  inviscid  and  the 
viscous  flows,  were  analyzed  and  compared  with  dual-wavelength  interferometric 
data  by  Liu  et  al  (Refs.  3,  4).  Satisfactory  agreement  was  obtained  for  the 
low-Mach-number  case,  Ms  = 13.1.  The  experimental  results  lie  between  the  frozen 
and  nonequilibrium  solutions.  However,  it  was  found  that  two-dimensional  effects 
are  significant  for  the  higher  Mach  nuzzfcer  case,  Ms  = 15.9*  The  phenomenon  of 
the  electron- cascade  front  moving  toward  the  wall  to  approach  the  translational 
shock  wave  was  found  more  prominent  for  stronger  shock  waves  by  Glass  et  al 
(Ref.  2).  The  reasons  for  this  premature  ionization  and  relaxation  process 
close  to  the  wall  are  still  unknown. 

The  present  sidewall  boundary-layer  results  were  obtained  by  using  the 
shock- structure  solutions  described  previously  as  the  edge  conditions  of  the 
boundary  layer.  The  difference  between  the  previous  model  by  Liu  et  al  (Ref.  4) 
without  taking  into  account  the  mutual  interactions  and  the  present  results  arises 
from  the  variations  of  freestream  flow  properties.  As  expected,  a change  in  the 
freestream  flow  conditions  will  cadse  changes  in  the  boundary-layer  flows. 

4.2  Comparison  of  Numerical  and  Experimental  Results  and  Discussions 

Figure  l4  shows  the  corresponding  first  order  and  present  results  for 
Ms  = 13.1,  Po  = 5»l6  torr  and  To  = 300  K together  with  the  experimental  data  for 
the  total-plasma-density  profiles,  p/p5,  where  6 denotes  the  conditions  outside 
the  boundary  layer.  Curve  A shows  the  first-order  results  for  the  nonequilibrium 
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boundary  layer;  curve  B shows  the  present  results  for  the  nonequilibrium  boundary 
layer  and  curve  C shows  the  first-order  results  for  the  frozen  boundary  layer. 

The  difference  between  the  first-order  and  the  present  coupled  results  is  mna.il 
since  it  was  shown  before  by  Liu  et  al  (Refs.  3,  *0  that  the  total  plasma- density 
profile  is  not  a sensitive  indicator  of  boundary-layer  structure.  However,  the 
difference  between  the  first-order  and  the  present  results  for  the  electron-number- 
density  profiles  can  be  readily  seen  in  Fig.  15,  even  though  the  difference  is 
small.  As  drawn  in  curve  B,  the  experimental  data  are  closer  to  the  present 
results.  This  shows  that  the  experimental  boundary-layer  profile  for  ne/neg  is 
closer  to  the  nonequilibrium  profile.  In  this  case  the  sidewall  boundary-layer 
profiles  are  measured  at  x = 9*5  cm  behind  the  shock  front.  The  relaxation  length 
is  about  9-°  cm  as  predicted  in  the  present  model  and  in  experiments.  Therefore, 
curve  B still  shows  strong  influence  from  the  electron- cascade  front  where  the 
profile  is  very  much  different  from  the  similar  profile  in  Liu  et  al  (Ref.  4) . 

The  same  tendency  is  observable  in  the  degree  of  ionization  profile  plotted  in 
Fig.  16.  It  is  evident  that  the  present  results  agree  better  with  the  experi- 
mental data  than  with  the  first-order  results  for  lower  shock  Mach  number,  where 
premature  ionization  and  recombination  close  to  the  wall  is  not  significant  in 
the  experiments  (see  Glass  et  al,  Ref.  2).  The  premature  ionization  is  discussed 
later  in  some  detail. 

Figure  17  shows  the  first-order  and  present  results  for  the  higher  shock 
Mach  nuuber  case  Ms  = 15*9  and  p0  = 5-10  torr  for  the  total  plasma  density  p/pg. 

The  agreement  between  the  present  result  and  experimental  data  is  not  as  good  as 
expected  in  the  lower  shock  Mach  number  case.  The  electron  number  density  profiled 
are  shown  in  Fig.  18.  The  present  results  (curve  B)  lie  between  the  nonequilibrium 
profile  ( curve  A)  and  the  frozen  profile  ( curve  C) , and  apart  from  the  experimental 
data.  The  experimental  data  differ  from  the  analytical  results  outside  the 
boundary  layer  edge  at  Y > 2 mm.  This  obviously  shows  the  two-dimensional  effect 
of  the  ionizing  shock-wave  structure  in  higher  shock  Mach  number  (see  Glass  et  al. 
Ref.  19). 

A tendency  similar  to  the  electron-number-density  profile  is  shown  in 
Fig.  19  for  the  degree  of  ionization  a/ag.  The  agreement  of  the  present  results 
with  the  experimental  data  is  worse . The  reason  for  the  discrepancy  between  the 
lower  and  higher  shock  Mach  nunber  cases  is  argued  as  follows.  It  was  shown  by 
Glass  et  al  (Ref.  2)  that  coupled  sinusoidal  disturbances  occur  at  higher  shock 
Mach  numbers  in  the  translational  shock  front  as  well  as  in  the  electron  cascade 
front  and  beyond.  Recently  the  flow  profiles  in  the  direction  normal  to  the 
wall  were  measured  by  Glass  et  al  (Ref.  19)  in  the  inviscid  flow  region.  It  was 
found  that  the  flow  is  nonuniform  in  the  inviscid  region  at  this  high  Mach  number 
Ms  = l6.  In  the  experiment,  the  electron- number  density  as  well  as  the  degree  of 
ionization  increased  with  distance  Y from  the  wall  and  did  not  reach  asymptotic 
values  at  the  outer  edge  of  the  boundary  layer  Y = l4  mm,  whereas  the  total-plasma 
density  decreased  from  the  wall  and  reached  its  asymptotic  value  at  the  edge  of 
the  boundary  layer.  However,  it  decreased  again  when  the  distance  from  the  wall 
was  greater  than  8 mm.  This  flow  nonuniformity  behind  a moving  shock  wave  at 
higher  shock  Mach  numbers  results  mainly  from  radiation  energy  transfer.  The 
radiation  energy  losses  are  prominent  in  the  inviscid  core  with  respect  to  the 
nonequilibrium  boundary-layer  flow  ( see  Fig . b6) . 

The  two-dimensional  effects  on  the  shock  structure  cannot  be  predicted 
simply  by  using  the  quasi-one-dimensional  models  as  used  in  the  previous  and 
present  analysis.  Further  analysis  such  as  using  the  unsteady  two-dimensional 
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model  may  be  tried  in  order  to  know  propagation  of  disturbances  induced  by  radiation- 
energy  transfer.  In  addition  to  this,  although  we  have  argued  about  the  optically 
thin  argon  plasma,  it  may  be  necessary  in  a future  study  to  take  into  account 
effects  of  reabsorption  of  the  radiation  energy  loss  in  the  inviscid  core  as  well 
as  in  the  boundary  layer.  This  effect  is  very  difficult  to  include  because  the 
reabsorption  coefficient  is  a f unction  of  the  complete  structure  of  the  radiation 
cooling  zone.  This  effect  was  not  included  in  the  present  analysis  (see  Liu  et  al, 
Ref.  3). 


Figures  20(a)  and  20(b)  show  experimental  contours  of  constant  degree 
of  ionization  a for  shocks  (Mg  = l6.5  and  Mg  = 13-6)  redrawn  from  Glass  et  al 
(Ref.  2) . The  solid  lines  indicate  the  outer  edge  of  the  degree  of  ionization 
boundary  layer  cc/a§  = 0.99,  obtained  from  the  present  sidewall  boundary-layer 
analysis  for  (a)  Mg  = 16.6,  pQ  = 4 .81  torr.  To  = 296  K and  ag  = 15. Tt>,  and  (b) 

Ms  = 13.1,  po  = 5.16  torr,  T0  = 300  K and  Og  = 6.4^.  These  small  discrepancies 
of  initial  conditions  between  analyses  and  experiments  do  not  lead  to  significant 
errors  in  comparing  both  results.  The  boundary- layer  thickness  for  degree  of 
ionization  reaches  a maximum  at  a distance  where  the  electron  number  density 
production  rate  He  is  a maximum,  and  then  temporarily  decreases.  The  reasons 
for  this  strange  behaviour  of  the  degree  of  ionization  can  be  explained  as  follows. 

Boundary- layer  flows  are  first  assumed  to  be  similar  just  behind  the 
shock  front.  Boundary-layer  flows  become  nonequilibrium  gradually  as  x increases. 
In  the  electron- cascade  front,  the  flow  properties  are  much  different  from  the 
similar  original  flow  properties  since  the  electron-nurber-density  production 
rate  and  the  radiation- energy  losses  are  greatest  there,  as  shown  in  Fig.  l6, 
for  example.  Consequently,  the  boundary- layer  profiles  are  also  changed  sharply 
with  distance  x near  the  electron- cascade  front.  Far  behind  the  electron- 
cascade  front,  the  boundary- layer  profiles  recover  to  the  weak  nonequilibrium 
profiles,  since  radiation-energy  losses  are  prominent  and  new  uniform  flow 
conditions  are  established,  for  example,  Fig.  4.  The  transition  of  the  boundary- 
layer  profile  from  a strong  nonequilibrium  profile  to  a new  similar  profile  occurs 
just  behind  the  electron  cascade  front.  In  addition  the  degree  of  ionization 
decreases  from  its  quasi- equilibrium  value  in  the  inviscid  core.  Therefore  the 
boundary-layer  thickness  of  the  degree  of  ionization  (a/ctg  = 0.99)  temporarily 
decreases,  as  shown  in  Fig.  20(a) . It  is  found  that  the  overshoot  in  the  a 
profile  across  the  boundary  layer  occurs  just  behind  the  electron- cascade  front 
for  the  following  reasons.  The  velocity  boundary  layer  as  well  as  the  boundary- 
layer  displacement  thickness  develop  almost  monotonically  with  distance  x by 
the  same  mechanism  as  in  the  flat-plate  boundary-layer  case  [see  Figs.  Bl(a)  and 
Bl(b)].  Their  growths  are  also  encouraged  by  acceleration  of  particle  velocity 
in  the  inviscid  core.  The  temperature  in  the  inviscid  core  decreases  so  sharply 
at  the  electron-cascade  front  that  the  temperature  boundary  layer  cannot  develop 
as  quickly  as  the  velocity  boundary  layer.  The  temperature  boundary  layer 
decreases  temporarily.  Consequently,  for  total  enthalpy  to  be  maintained 
across  the  boundary  layer,  the  temperature  profile  across  the  boundary  layer 
should  have  an  overshoot.  As  a result,  this  temperature  ever shoot  across  the 
boundary  layer  induces  the  overshoot  of  degree  of  ionization  profile  at  the 
electron- cascade  front. 

It  is  clear  that  this  overshoot  significantly  affects  the  boundary- 
layer  thickness  of  the  degree  of  ionization,  63,  defined  by  Eq.  (28d).  In 
Fig.  B5,  variations  of  Gq  with  distance  x are  shown  for  shocks  with  Mg  = 15*9 
and  13.I.  It  is  worth  noting  that  due  to  this  overshoot  profile,  §3  can  be 
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a maximum  near  the  electron- cascade  front  in  the  case  of  the  nonequilibrium 
boundary  layer.  However,  in  the  frozen  boundary  layer,  the  overshoot  of  a is 
not  found.  For  details  see  Appendix  B. 

Arrows  with  figures  along  a/otg  = 0.99  curves  in  Fig.  20(a)  show  the 
degree  of  ionization  (%)  at  the  marked  points.  The  value  of  a along  the  curve 
a/ofg  = 0-99  increases  monotonically  from  0$  at  the  shock  front  to  4%  at  x = 1.25 
cm.  Then  it  suddenly  increases  across  the  electron- cascade  front  and  reaches  a 
maximum  of  a = 15.7%.  It  is  shown  that  experimental  a contours  indicate  the  a 
overshoot  at  the  electron- cascade  front.  The  agreement  of  a between  the  present 
analysis  and  experiments  is  very  good  at  and  behind  the  electron- cascade  front; 
but  it  is  poor  in  the  incubation  region  between  the  shock  front  and  the  electron- 
cascade  front,  since  the  2%  contour  intersects  with  the  a/ 0£g  = 0.99  curve  where 
the  predicted  values  are  less  than  1%.  The  reason  for  this  discrepancy  is  not 
clear.  The  same  remarks  can  be  applied  to  the  lower  shock  Mach  number  case, 

Ms  - 13.1  as  drawn  in  Fig.  20(b). 

Agreement  between  experimentally  obtained  a contours  and  predicted 
a/ag  = 0.99  curves  is  very  good.  It  is  worth  noting  that  the  premature 
ionization  near  the  sidewall  is  merged  with  the  predicted  degree  of  ionization 
boundary  layer. 

It  is  known  that  far  behind  the  electron-cascade  front  the  degree  of 
ionization  and  electron-number  density  in  experiments  do  not  approach  their 
asymptotic  values  but  still  Increase  beyond  the  outer  edge  of  the  boundary 
layer  for  higher  shock  Mach  number  Ms  = 15.9.  Glass  et  al  (Ref.  17)  studied 
these  phenomena  and  concluded  that  a two-dimensional  effect  is  mainly  due  to 
the  radiation  energy  transfer  mechanism.  The  details  are  still  far  from  clear. 
Experimental  indications,  for  example  in  Fig.  20(a),  show  that  the  degree  of 
ionization  profile  across  the  boundary  layer  at  the  electron- cascade  front 
x <=*  17  mm  initially  increases  with  distance  Y from  the  sidewall.  Then  it 
reaches  the  quasi-equilibrium  peak  value  Oft  = 15.8%.  Finally  it  gradually 
decreases  in  the  inviscid  core  Y > 2 mm  since  the  experimental  a contours  are 
concave  to  the  shock  front.  It  is  found  that  this  profile  shows  a completely 
different  tendency  from  that  far  behind  the  electron-cascade  front  at  x = 18  cm 
as  shown  by  Glass  et  al  (Ref.  19)  or  in  Fig.  19.  The  reasons  for  this  differ- 
ence are  explained  as  follows.  The  a profile  across  the  shock  tube  is  affected 
by  the  boundary- layer  growth,  radiation  energy  losses  and  ionization  and 
recombination  rate  processes.  Since  these  effects  themselves  and  their  coupled 
effects  vary  significantly  with  distance  x,  it  is  quite  acceptable  to  assume 
that  the  experimental  Q£  profile  measured  at  the  electron- cascade  front  is  very 
different  from  that  measured  far  behind  the  electron- cascade  front. 

In  an  actual  unsteady  two-dimensional  flow,  especially  behind  the 
electron- cascade  front,  disturbances  generated  by  the  boundary- layer  growth, 
radiation  energy  losses  and  ionization  and  recombination  processes  are  non- 
uniformly  propagating  along  the  characteristic  surfaces.  These  nonuniformly 
propagating  disturbances  interact  with  the  electron- cascade  front  as  well  as 
the  shock  front.  However,  it  is  far  from  clear  how  these  disturbances  interact 
two- dimensionally  with  the  electron- cascade  front  and  how  they  are  related  to 
the  sinusoidal  configuration  of  a perturbed  shock  front,  as  found  by  Glass  et  al 
(Ref.  2).  Although  the  premature  ionization  near  the  sidewall  is  found  to  be 
merged  with  the  degree  of  ionization  boundary  layer,  their  detailed  interpretations 
are  still  far  from  being  completed.  It  will  be  an  important  future  research  task 
to  understand  their  origin  and  interactions  with  ionizing  argon  shock-tube  flows. 
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Deckker  et  al  (Ref.  20),  applying  the  hot  wire  anemometry  to  shock- 
tube  flows,  showed  that  flow  property  distributions  are  not  uniform  close  to 
the  shock  front  even  at  very  low  shock  Mach  number  Ms  - 1.16  (the  pressure 
ratio  across  the  shock  front  Ppi  = 1.57) • Velocity  and  temperature  profiles 
are  found  to  be  different  from  those  profiles  predicted  by  the  similar  boundary- 
layer  analyses.  Although  a critical  examination  of  the  reproducibility  of  their 
shock- tube  experiments  and  disturbing  effects  of  the  probe  on  the  subsonic 
shock-tube  flow  is  needed,  the  question  arises  as  to  which  mechanisms  cause 
the  flow  nonuniformity  close  to  the  shock  front.  They  concluded  that  the 
flow  nonuniformity  may  be  due  to  the  unsteady  effects  of  the  shock-tube  flow, 
especially  due  to  the  propagation  of  disturbances  generated  at  the  contact 
front. 

Further  experiments  and  more  sophisticated  analyses  applying  the 
unsteady  two  or  three-dimensional  Navi er- Stokes  equations  to  the  ionizing 
argon  shock- tube  flow  are  necessary  to  understand  the  flow  nonuniformity. 


5.  QUASI-STEADY  FLAT-PLATE  BOUNDARY-LAYER  FLOWS 
5 .1  General  Discussions 


The  first-order  results  of  quasi- steady  flat-plate  boundary-layer 
flows  generated  by  a strong  shock  wave  in  a shock  tube  for  ionizing  argon 
were  given  and  compared  with  dual -wavelength  interferometric  data  by  Liu  et  al 
(Ref.  3).  Satisfactory  agreement  was  obtained  between  analysis  and  experi- 
mental total  plasma-density  and  electron- number  density  profiles  for  shocks 
Mg  = 13.1  and  15-9* 

Two  cases  reported  previously,  Ms  = 12.8,  p0  = 5*01  torr  and  To  = 

297  K,  and  Ms  = l6.6,  po  = 4.8l  torr  and  T0  = 298  K,  are  re-analysed  by  using 
the  present  method.  The  difference  between  the  previous  result  by  Liu  et  al 
(Ref.  3)  and  the  present  analysis  arises  from  the  changes  in  the  free- stream 
flow  conditions  which  result  from  mutual  interactions  between  inviscid  and 
viscous  flows. 

5.2  Conparison  of  Numerical  and  Experimental  Results  and  Discussions 

Figure  21  shows  comparisons  between  analyses  and  experimental  data 
for  the  case  with  Ms  = l6.6,  pQ  = 4.8l  torr  and  T0  = 297  K for  the  total- 
plasma-densi L j . The  difference  between  the  present  analysis  and  the  first- 
order  results  is  quite  small.  The  corresponding  results  for  the  electron- 
number  density  ne/n^  and  degree  of  ionization  profiles  Of/ag  are  plotted  in 
Figs.  22  and  23,  respectively.  A bump  appears  in  the  experimental  data  for 
the  electron-number  density  profile  and  in  the  present  analysis,  but  the 
maximum  deviation  from  the  asymptotic  value  of  ne/ne6  -»1  is  very  small,  at 
most  2 The  disagreement  between  analysis  and  experiments  for  the  electron- 
number  density  results  mainly  from  the  fact  that  the  boundary-layer  flow  is 
assumed  to  be  quasi-steady  while  in  the  experiments  it  is  unsteady  due  to 
radiation-energy  losses.  That  is,  the  electron-number  density  and  degree 
of  ionization  are  slowly  decreasing  with  distance  from  the  electron-cascade 
front  and  the  pressure,  particle  velocity  and  total  plasma  density  are  slightly 
increasing  with  distance  from  the  electron-cascade  front,  as  discussed  in 
Chapter  3.  However,  variations  of  these  flow  properties  with  distance  x were 
not  considered  in  the  present  flat-plate  boundary- layer  flow  in  order  to  satisfy 
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the  steady- state  conditions.  The  boundary-layer  displacement  thickness  growths 
are  shown  with  distance  x in  Fig  2k  for  both  nonequilibrium  and  frozen  boundary- 
layer  flows.  As  expected  in  the  steady  flat-plate  flow,  the  displacement  thick- 
ness grows  almost  proportional  to  the  square  root  of  the  distance  from  the 
leading  edge  of  the  flat-plate,  whereas  in  the  sidewall  boundary  layer,  bumps 
appear  at  the  electron- cascade  front.  The  boundary-layer  displacement  thickness 
in  the  frozen  flow  is  thinner  than  that  in  the  nonequilibrium  flow,  although 
the  Reynolds  numbers  referring  to  the  inviscid  core  flow  properties  are  the  same. 
The  reason  is  that  in  the  nonequilibrium  boundary-layer  flow,  kinetic  energy 
is  spent  for  internal  electronic  excitation  of  argon  atoms  as  well  as  for 
ionization  processes.  As  a result,  a temperature  gradient  and  other  flow 
property  gradients  across  the  boundary  layer  are  comparatively  smoothed  out. 

For  example,  Figs.  23  and  2k  show  the  electron  number  density  and  degree  of 
ionization  profiles,  respectively.  Therefore  the  nonequilibrium  boundary- layer 
displacement  thickness  is  thicker  than  that  of  the  frozen  flow. 

In  Table  2 comparisons  of  the  uniform  conditions  between  the  first- 
order  results  and  the  present  results  are  tabulated.  The  differences  between 
these  results  are  not  significant.  The  local  flow  Mach  numbers  M§  are  increased 
by  about  3 % for  Ms  = l6.6  and  by  about  5$  for  Ms  = 12.8.  For  the  latter  case, 

Ms  = 12.8,  the  nonequilibrium  and  frozen  flow  profiles  for  total  plasma  density 
p/pg  and  for  electron-number  density  are  shown  in  Figs.  25  and  26,  respectively, 
together  with  experimental  data.  A small  improvement  is  found.  The  experimental 
results  for  the  electron-number  density  lie  between  the  analytical  nonequilibrium 
and  the  frozen- flow  profiles.  The  corresponding  results  for  the  degree  of  ioniza- 
tion are  shown  in  Fig.  27. 

Figure  28  shows  variations  of  nonequilibrium  and  frozen  boundary- 
layer  displacement  thicknesses  with  distance  x for  lower  shock  Mach  number 
Ms  = 12.8.  It  is  found  that  the  frozen  boundary- layer  displacement  thickness 
is  mauch  thinner  than  that  of  the  nonequilibrium  flow.  Effects  of  ioniza- 
tion and  recombination  processes  are  more  significant  in  the  lower  shock  Mach 
number,  as  shown  in  Figs.  26  and  27-  Comparing  the  result  for  Ms  = 12.8  with 
the  higher  shock  Mach  number  Ms  = l6.6,  it  is  found  that  the  growth  of  the 
boundary- layer  displacement  thickness  for  the  lower  shock  Mach  number  is  much 
slower  than  that  for  the  higher  shock  Mach  number.  The  Reynolds  numbers 
referring  to  the  inviscid  freestream  properties  are  different  in  each  case. 

It  is  shown  that  effects  of  the  mutual  interactions  between  the 
inviscid  and  viscous  flows  on  the  quasi- steady  flat-plate  boundary- layer  flows 
are  small.  However,  the  flow  model  used  in  the  present  analysis  is  a quasi- 
steady one.  Unsteady  effects  appearing  in  the  experiments  are  not  considered 
in  the  analysis.  Therefore,  in  order  to  obtain  better  agreement  between  analysis 
and  experiments,  the  unsteady  two-dimensional  Navier-Stokes  equations  should  be 
solved  with  proper  unsteady  inviscid  freestream  conditions.  They  are,  of  course, 
one  of  the  important  research  tasks  for  the  future. 


6.  DISCUSSIONS  AND  CONCLUSIONS 

An  analysis  is  presented  of  the  mutual-interaction  effects  of  ionizing- 
argon  shock-wave  structure  and  shock-tube  sidewall  as  well  as  flat-plate  boundary- 
layer  flows.  The  effective  quasi-one-dimensional  shock-wave- structure  equations 
for  a tube  of  finite  area  are  derived  from  the  flow  area-averaged  concept.  Flow 
nonuniformities  resulting  from  the  boundary-layer-displacement  thickness  and 
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other  defined  thicknesses,  as  well  as  the  wall  dissipation  terms,  are  included 
in  the  shock- wave- structure  analysis.  The  results  are  compared  with  a previous 
first-order  analysis  and  interferometric  data  at  nominal  shock  Mach  n ushers  of 
13  and  16  at  an  initial  pressure  of  5 torr.  It  is  shown  that  a perfect-gas 
analysis  based  only  on  displacement  thickness  overestimates  the  reduction  in 
shock- wave- structure  thickness.  However,  a real-gas  analysis  shows  that  the 
actual  shock-wave  structure  thickness  is  closer  to  the  idealized  one- dimensional 
case.  In  addition,  the  interferometric  measurements  of  the  electron-number 
density  and  total-plasma  density  shock- wave- structure  profiles  as  well  as  the 
derived  degree  of  ionization  profiles  are  in  much  better  agreement  with  the 
present  analysis  than  those  done  previously.  The  same  remarks  apply  to  the 
ionizing  sidewall  as  well  as  flat-plate  boundary- layer  interferometric  profiles. 
Agreement  of  the  experimental  data  with  the  present  analysis  is  very  good  for 
both  the  shock-wave  structure  and  the  boundary  layer  at  Mg  ~ 13  and  good  at 
Mg  ~ 16.  It  is  found  that  the  premature  ionization  near  the  sidewall  is  merged 
with  the  degree  of  ionization  boundary- layer  thickness. 

Further  investigation  based  on  the  time -dependent  Navier-Stokes 
equation  may  be  required  in  order  to  predict  the  two-dimensional  effects  on 
shock-wave  structure  and  the  unsteady  effects  on  induced  flat-plate  boundary- 
layer  flows. 
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Table  1(a) 


uasi-equilibrium  peak  values  of  a,  n and  p and  relaxation  lengths 


The  first  iteration  results  at  M_  = 13* ** ***1>  P„  = 5-16  torr  and  T^  = 3 


Mg  " 13.1 


First  order 


Non- 

equilibrium 


Frozen*** 


1 D ! 

*E 

cm 

cm 

00 

10.00 

6.43 

6.63 

12.86 

7.93 

8.93 

25.71 

8.88 

6.43 

7.88 

12.86 

8.98 

25.71 

9.48 

0.0590 


0.0711 


* Results  by  higher  iteration  number 

**  Nonequilibrium  boundary- layer  profile 

***  Frozen  boundary- layer  profile 


xl017/cm3 


0.566 


.771 


Pe 

-4  . 3 
xlO  g/  can 


0.650 


0.719 

0.680 
0.660  * 

O.661 


0.681 

0.660 

0.650 


Table  l(b) 


Quasi-equilibrium  peak  values  of  a,  ng  and  p and  relaxation  lengths. 
The  first  iteration  results  at  Mg  = 15 -9 > PQ  = 5*10  torr  and  TQ  = 297  K. 


Ms  =15-9 

1 

D 

*E 

“E 

neE 

PE 

can 

cm 

xlOl8/cm3 

-4  , 3 
xlO  g/  cm0 

First  order 

00 

2.15 

0.1361 

0.168 

0.821 

6.43 

1.83 

0.1414 

0.182 

0.852 

Non- 

equilibrium 

** 

12.86 

2.03 

2.03 

0.1363 

0.1359 

0.170 

0.170 

0.828 
0.830  * 

25.71 

2.12 

0.1355 

0.168 

0.824 

6.43 

1.93 

0.1424 

0.182 

0.848 

Frozen*** 

12.86 

2.03 

0.1371 

O.D3 

0.825 

25.71 

2.13 

0.1362 

0.170 

0.823 

* Results  by  higher  iteration  number 

**  Nonequilibrium  boundary-layer  profile 

***  Frozen  boundary- layer  profile 


Table  2 


Initial  conditions  and  freestream  conditions  for  flat-plate  boundary- layei’  flow. 
Con^arison  of  present  conditions  with  those  in  Ref.  3- 


Ms 

Po  (torr) 

*„  (’« 

l6.6 

4.81 

296 

12.8 

5.01 

297 



Ref.  3 

Present 

Ref.  3 

Present 

Ug  (cm/s) 

4 .86xl05 

4 .97xl05 

3-53X105 

3-59X105 

P&  (torr) 

2025 

2160 

1200 

1189 

Ta6  <“*> 

l . 049x10^ 

1. 043x10** 
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FIG.  1(a)  COMPARISON  OF  BOUNDARY-LAYER  DISPLACEMENT-THICKNESS  GROWTHS  FOR  THREE  ITERATIONS 
FOR  Ms  = 13.1,  pQ  « 5.16  TORR  AND  TQ  = 300  K. 
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FIG.  1(b)  COMPARISON  BOUNOARY-LAYER  DISPLACEMENT-THICKNESS  GROWTHS  FOR  TWO  ITERATIONS  FOR 

M * 15.9,  p <=  5.10  TORR  ANO  T » 398  K. 
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FIG.  2 VARIATIONS  OF  NORMALIZED  ELECTRON-NUMBER  DENSITY  ng/neE  WITH  OlSTAflCE  X FOR  M$  * 13. 1, 
pQ  = 5.16  TORR  AND  TQ  = 3C0  K.  SUBSCRIPT  E OENOTES  QIIAS I -EQUI L IBR I UM  VALUE.  CURVE 
A:  FIRST-ORDER  RESULTS  WITHOUT  MUTUAL  INTERACTIONS;  CURVE  B:  ENOMOTO'S  MODEL;  CURVE 
C:  FIRST-ITERATIVE  RESULT;  CURVE  D:  FINAL  ITERATIVE  RESULT.  ARROWS  SHOW  TERMINATION 
OF  RELAXATION. 


FIG.  3 VARIATIONS  OF  NORMALIZED  TOTAL-PLASMA-DENSITY  PROFILES  p/c>c  WITH  DISTANCE  X FOR 
MS  ' ,3-1,  Po  ' 5-16  T0RR  AND  To  ‘ 300  K>  (f0R  0TMCR  NOTATIONS,  SEC  FIG.  2.) 
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FIG.  4 VARIATIONS  OF  NORMALIZED  DEGREE  OF  IONIZATION  PROFILES  a/*c  WITH  DISTANCE  X FUR 
Ms  - 13.1.  pQ  » 5.16  TORR  AND  T0  • 300  K.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  2.) 


FIG.  5 VARIATIONS  OF  NORMALIZED  PRESSURE  P/P0  WITH  DISTANCE  X FOR  M,.  ■ 13.1,  P0  • 5.16  TORR 
ANU  T • 300  K.  CURVE  A:  FIRST-ORDER  RESULTS  WITHOUT  MUTUAL  INTERACTIONS;  CURVE  B; 
ENOMOTO'S  MODEL;  CURVE  0:  FINAL  ITERATIVE  RESULTS  OF  PRESENT  MOOEl. 
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FIG.  8 VARIATIONS  OF  NORMALIZED  ELECTRON-NUMBER -DENSITY  PROFILES  ne/ne£  WITH  DISTANCE  X FOR 
M • 15.9,  p * 5.1  TORR  AND  T ■ 298  X.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  2.) 
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FIG.  9 VARIATIONS  OF  NORMALIZED  TOTAL-PLASMA-DENSITY  PROFILES  p/»E  WITH  DISTANCE  X FOR 
Ms  ’ I5-9,  po  ’ SJ  T0RR  W<0  To  ' 298  *•  (F0R  0THER  r,0TATI0NS-  SEE  FIG.  2.) 
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FIG.  11  VARIATIONS  OF  NORMALIZED  PRESSURE  P/PQ  WITH  DISTANCE  * FOR  Ms  * 15.9,  pQ  * 5.10  TORR 
AND  T0  * 298  K.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  5.) 
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FIG.  13  VARIATIONS  OF  NORMALIZED  ELECTRON-NUHBER-DEIISITY  PROFILES  WITH  DISTANCE  X 

FOR  VARIOUS  SHOCK-TUBE  HYDRAULIC-DIAMETER  D FOR  CASE  KITH  M$  ■ 15.9,  P0  * 5.1  TORR 
AND  Tq  - 298  K.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  7.) 
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FIG.  14  COMPARISON  OF  ANALYTICAL  AMO  EXPERIMENTAL  NORMALIZED  PROFILES  OF  PLASMA  DENSITY 
P/Pt  WITH  DISTANCE  Y IN  THE  SIDEWALL  BOUNDARY  LAYER  IN  ARGON  AT  X • 9.5  CM, 

Ms  - 13.1,  PQ  * 5.16  TORR  AND  TQ  • 300  K.  CURVE  A:  FIRST-ORDER  RESULT;  CURVE  B; 
PRESENT  RESULT;  CURVE  C:  FROZEN  SOLUTION. 


FIG.  15  COMPARISON  OF  ANALYTICAL  AIIO  EXPERIMENTAL  NORMALIZED  ELECTRON-NUttER-DENSITY  PROFILES 

n9/nft  WITH  DISTANCE  Y IN  THE  SIDEWALL  BOUNDARY  LAYER  IN  ARGON  AT  X ■ 9.5  CM,  H$  • 13.1, 

p*  -*5.16  TORR  AND  T„  • 300  K.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  14.) 
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FIG.  18  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  PROFILES  OF  ELECTRON  HUMBER  DENSITY 
Vne«  W!TH  0ts™ce  r ,N  TH£  SWCHALL  BOUNDARY  LAYER  AT  X • 18  CM,  M$  - 15.9,  pQ  * 5. 1 
TORR  ANU  TQ  ■ 298  K.  (FOR  OTHER  NOTATIONS,  SEE  TIG.  14.) 
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FIG.  19  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  PROFILES  OF  DEGREE  OF  IONIZATION 
a/os  WITH  DISTANCE  Y IN  THE  SIDEWALL  BOUNDARY  LAYER  AT  X • 18  CM,  Mt  - 15.9,  P0  - 5.1 
TORR  AND  T0  • 298  X.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  14.) 


FIG.  20(a)  CONTOURS  OF  CONSTANT  DEGREE  OF  IONIZATION  IN  RELAXATION 

REGION  CLOSE  TO  WALL  IN  ARGON  AT  M$  = 16.5,  pQ  • 5.12' TORR 
AND  TQ  - 297  K (SEE  REF.  2).  SOLID  LINE:  CONTOUR  OF 
DEGREE  OF  IONIZATION  BOUNDARY  LAYER  a/a&  = 0.99  AT 

Ms  * ™.6,  P0  * 4-81  T0RR  AND  To  * 298  K- 


FIG.  20(b)  CONTOURS  OF  CONSTANT  DEGREE  OF  IONIZATION  IN  RELAXATION 

REGION  CLOSE  TO  WALL  IN  ARGON  AT  Hs  * 13.6,  pQ  ■ 5.09 

AND  TQ  - 296.7  K (SEE  REF.  2).  SOLID  LINE:  CONTOUR 

OF  DEGREE  OF  IONIZATION  BOUNDARY  LAYER  <x/a6  AT  M$  - 13.1, 

- 5.16  TORR  AND  T - 300  K. 
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FIG.  21  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  PROFILES  OF  PLASMA  DENSITY  p/o{ 
WITH  DISTANCE  Y FOR  FLAT-PLATE  BOUNDARY  LAYER  AT  X = 14  CM,  M$  = 16.6,  PQ  * 4.81  TORR 


AND  T = 296  K.  CURVE  A:  PRESENT  NONEQUILIBRIUM;  CURVE  B:  PRESENT  FROZEN;  CURVE  C; 
0 

FIRST-ORDER  NOHEQUILIBRIUM;  CURVE  D:  FIRST-ORDER  FROZEN. 


FIG.  22  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  ELECTRON-NUMBER-DENSITY  PROFILE 

n /n  WI™  DISTANCC  y F0R  fLAT-PLATE  BOUNDARY  LAYER.  (INITIAL  CONDITIONS  AND  NOTATIONS 
ARE  GIVEN  IN  FIG.  21.) 
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FIG.  23  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMAL I2ED  PROFILE  FOR  DEGREE  OF  IONIZATION 

a/a,  WITH  DISTANCE  Y FOR  FLAT-PLATE  BOUNDARY  LAYER.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  21.) 
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FIG.  25  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  PROFILE  FOR  PLASMA  DENSITY 
o/p{  FOR  FLAT-PLATE  BOUNDARY  LAYER  AT  X « 14  CM,  M$  * 12.8,  pQ  ' 5.01  TORR  AND 
T0  - 297  K.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  21.) 


FIG.  26  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  PROFILE  FOR  ELECTRON  NUMBER 
DENSITY  ne/n  !’  THE  FLAT-PLATE  BOUNDARY  LAYER  AT  X * 14  CM,  = 12.8,  p0  ■ 5.0) 

TORR  AND  Tq  X.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  21.) 


FIG.  27  COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  NORMALIZED  PROFILE  FOR  DEGREE  OF  IONIZATION 
a./ c{  FOR  FLAT-PLATE  BOUNDARY  LAYER  AT  X = 14  CM,  M$  = 12.8,  pQ  = 5.01  TORR  AND  Tq  * 
297  K.  (FOR  OTHER  NOTATIONS,  SEE  FIG.  21.) 


FIG.  28  VARIATIONS  OF  FLAT-PLATE  NONEQUILIBRIUM  AND  FROZEN  BOUNDARY-LAYER  DISPLACEMENT- 
THICKNESSES  WITH  DISTANCE  X FOR  M - 12.8,  p ' 5.01  TORR  AND  T - 297  K 
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APPENDIX  A 


AREA  AVERAGED  QUASI-ONE-DIMENSIONAL  FLOW  EQUATIONS 


In  the  shock- fixed  coordinates,  the  area-averaged  mass  conservation 
equation  (Eq.  18)  is  rewritten  as  follows.  Let  the  operator  < > be  introduced 
in  an  area- averaging  process,  and  let  A and  L be  the  cross-sectional  area  and 
perimeter  of  a shock  tube,  respectively.  A schematic  profile  is  shown  in  Fig. 
Al. 


<pu> 


pudA 


i r 
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= j jpu  J dA  + L J (p’u*  - pu)dyj- 


(Al) 


b . layer 


i pu(A  - L6») 


where  the  primed  quantities  represent  variables  across  the  boundary  layer  and 
6*  is  the  boundary-layer-displacement  thickness  defined  by 


6*  = 


'u' 


1 - V ,dy 


(A2) 


b .layer 

which  is  negative,  as  expected,  p and  u are  the  uniform  values. 

Here  inside  the  boundary  layer,  the  mass  conservation  equation 


(A3) 


holds.  If  we  integrate  Eq.  (A3)  across  the  boundary  layer,  we  have  the  following 
relation: 


(p'v'h  = -&f  p'u' 


dy 


(A4) 


b .layer 

where  5 denotes  the  condition  at  the  outer  edge  of  the  boundary  layer. 


For  an  arbitrary  quantity  Q,  its  convective  derivative  averaged  across 
the  shock  tube  cross-sectional  area  is  given,  taking  into  account  Eqs.  (A^)  and 
(A3),  by 


A-l 


-*’•  - I 


s <PUQ>  + <PVQ> 

= J PUQ(A  - &*D  + | ^ f (P'u Q'  - p'u'Q)dy 

b .layer 

= I puQ^A  ~ + 6Q^  (A5) 

where  &'  is  defined  by 

vj  i9:(1-t:)ay  (*> 

b . layer 

This  is  the  commonly  used  boundary- layer  thickness  for  Q,  which  represents 
u,  a,  H and  OfTe  corresponding  to  the  momentum,  electron  mass,  energy  and 
electron  energy  conservation  equations,  respectively. 

Here  we  rewrite  the  boundary-layer  thickness  for  Q as 


6Q  = 5*  + 6^  (A7) 

-I  0-^)* 

b . layer 

Finally  we  have  an  area-averaged  quasi-one-dimensional  expres-  - 
sion  as  follows: 


<puQ>  = i puQ(A  - 5qL) 


(a8) 
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APPENDIX  B 


BOUNDARY  LAYER  THICKNESSES  AND  WALL  DISSIPATION  TERMS 


Once  various  boundary- layer  thicknesses  and  wall  dissipation  terms 
are  given,  we  can  determine  the  shock-wave  structure  in  ionizing  argon  by 
making  use  of  these  mutual  interaction  terms.  The  boundary-layer  thicknesses 
as  well  as  wall  dissipation  terms  used  in  the  present  analyses  are  presented 
here. 


In  Figs.  Bl( a)  and  Bl(b) , variations  of  the  boundary- layer  displace- 
ment thickness  -5*  and  the  velocity  thickness  -&u  with  distance  x are  shown 
for  shocks  with  shock  Mach  numbers  Ms  = 13.1  and  Ms  = 15 -9>  respectively.  The 
thicknesses  are  defined  by 


6* 
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p'u' 

pu 


dy  < 0 
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ul 
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dy  < 0 


(Bl) 


where  primed  quantities  show  those  across  the  boundary  layer.  5*  and  &u  are 
negative  in  the  shock- fixed  coordinate,  as  expected.  It  is  noted  that  the  bump 
appears  in  both  nonequilibrium  and  frozen  boundary-layer  flows  at  the  electron- 
cascade  front  where  the  electron  number  density  production  rate  ne  is  greatest. 
Therefore  these  boundary-layer  thickness  growths  are  no  longer  monatonically 
increasing  as  seen  in  the  similar  boundary  layer  analysis  for  perfect  argon. 

The  velocity  boundary- layer  thicknesses  -&u  temporarily  decrease  at  the 
electron- cascade  front  due  to  the  sharp  changes  of  boundary-layer  profiles. 

The  phenomena  are  even  greater  in  higher  shock- wave  Mach  numbers. 

Figure  B2  shows  the  variation  of  -6*/«^c  with  distance  x for  shock- 
waves  with  Ms  = 13.1  and  15-9.  In  the  similar  laminar  boundary- layer  analysis 
for  perfect  argon,  the  sidewall  boundary- layer  displacement  thickness  -6*  is 
growing  with  «7x.  Figure  B2  shows  roughly  a coefficient  of  the  displacement 
thickness  growth,  which  is  proportional  to  the  square  root  of  the  Reynolds 
number  referring  to  the  uniform  flow  condition.  The  value  of  -5*/s/x  is 
initially  constant  in  the  incubation  region  between  the  shock  front  and  the 
electron- cascade  front  and  then  suddenly  changes  to  another  quasi-equilibrium 
value  across  the  electron- cascade  front.  The  reasons  for  this  are  as  follows. 
Pressure  changes  roughly  from  a constant  value  in  the  incubation  region  to 
another  constant  value  behind  the  electron- cascade  front.  In  other  words, 
the  flow  Reynolds  number  changes  suddenly  across  the  electron- cascade  front. 
Values  of  -6 at  x = 0 (-&*/*Vx  = 0.1006  at  Ms  = 13*1  and  -d*/\/x  = 0.1064 
at  Ms  = 15-9)  are  important  because  these  values  are  used  in  calculating  shock- 
wave  structures  by  Enomoto's  model.  This  model  considers  only  the  boundary- 
layer  displacement  thickness  which  is  obtained  by  a similar  boundary-layer 
analysis.  In  Fig.  B3,  variations  of  momentum  thickness  &m  with  distance  x 
are  shown  for  shocks  Ms  = 13 •!  and  15 *9 > 


B-l 


C1  - 


(B2) 


Solid  lines  and  dashed  lines  show  nonequilibrium  and  frozen  boundary  layer, 
respectively.  Inflection  points  appear  at  the  electron- cascade  front  owing 
to  the  same  reasons  as  explained  in  Figs.  Bl(a)  and  Bl(b) . Differences 
between  nonequilibrium  and  frozen  boundary- layer  flows  are  found  to  be  very 
small. 

Figure  shows  the  energy  dissipation  thickness  6e  for  shocks 
Ms  = 13.1  and  15-9>  where 

6e  = 6**J'£5r(1-7-)4ir  (B3> 

Curve  A ( — -)  and  curve  B ( — ) correspond  to  nonequilibrium  and  frozen 
boundary- layer  flows,  respectively.  The  same  remarks  made  for  Fig.  B3  can  be 
applied  to  this  case. 

Figure  B5  shows  variations  of  the  degree  of  ionization  thickness 
6ct  with  distance  x for  shocks  Ms  = 13*1  and  15 -9 > where 

6a  = s*  +f^r  C1  ■ 5r)  ^ (B4) 

Curve  A ( ) and  curve  B ( — ) correspond  to  nonequilibrium  and  frozen 

boundary-layer  flows,  respectively.  As  mentioned  previously,  overshoots 
of  the  a profile  occur  across  the  nonequilibrium  boundary  layer  just  behind 
the  electron- cascade  front  as  shown  in  Figs.  l6  and  19 . They  do  not  appear 
in  the  frozen  boundary-layer  flow.  Consequently,  ^2  in  the  nonequilibrium 
boundary-layer  flow  can  be  positive  around  the  electron- cascade  front  and  is 
a mavimi-m  at  a distance  where  the  electron  number  density  production  rate  is 
greatest.  On  the  other  hand,  E«  in  the  frozen  boundary-layer  flow  is  always 
negative  and  monotonically  decreases  with  distance  x.  It  is  worth  noting 
that  a/ag  = 0.99  contours  in  Figs.  20(a)  and  20(b)  are  similar. to  the  results 
obtained  for  the  degree  of  ionization  boundary-layer  thickness  in  the  nonequi- 
librium boundary-layer,  which  is  defined  as 

/*£  (>-£)« 

Figure  B6  shows  variations  of  the  radiation  energy  loss  thickness 
6r  with  distance  x for  shocks  Ms  = 13*1  and  15 *9 > where 

wC1-!)'*  <b5) 

Curves  A ( ) and  B(  — ) correspond  to  nonequilibrium  and  frozen  boundary- 

layer  flows,  respectively.  5r  in  the  nonequilibrium  boundary  layer  is 


B-2 


always  positive.  This  fact  indicates  that  radiation  energy  losses  inside 
the  nonequilibrium  boundary  layer  are  smaller  than  those  in  the  inviscid 
core.  6r  is  a maximum  at  the  electron-cascade  front  where  the  averaged 
radiation  energy  loss  inside  the  nonequilibrium  boundary  layer  is  smallest. 

The  radiation  energy  losses  in  the  frozen  boundary  layer  are  larger  than 
those  in  the  inviscid  core,  since  5r  is  negative,  as  shown  in  Fig.  B6.  In 
the  nonequilibrium  boundary  layer,  kinetic  energy  is  transformed  into  excita- 
tion of  the  internal  electronic  states,  and  ionization  rate  processes.  Hence 
the  temperature  decreases  more  than  in  the  frozen  boundary  layer  where  the 
ionization  and  recombination  processes  are  assumed  completely  frozen.  There- 
fore, the  radiation  energy  losses  are  more  prominent  in  the  frozen  boundary- 
layer  flow  than  in  the  nonequilibrium  boundary-layer  flow.  However,  throughout 
the  previous  and  present  analysis,  the  effect  of  re-absorption  of  radiation 
energy  in  the  radiation  cooling  zone  is  not  considered. 

Figure  B7  shows  variations  of  the  nonequilibrium  electron- number 
density  thickness  &ne  with  distance  x for  shocks  Ms  = 13*1  and  15. 9>  where 


As  shown,  there  is  a maximum  at  the  electron- cascade  front.  The 
electron  number  density  in  the  inviscid  core  is  a maximum.  The  electron 
number  density  profile  across  the  nonequilibrium  boundary  layer  near  the 
electron- cascade  front  is  shown  in  Fig.  15  where  the  profiles  are  compared 
at  x = 9.5  cm  behind  the  shock  front  and  the  experimental  relaxation  length 
is  about  9.0  cm.  In  Fig.  l8  similar  comparisons  are  made  at  x = 18.0  cm 
behind  the  shock  front  and  the  experimental  relaxation  length  is  2.0  cm. 
There  is  evidence,  from  Figs.  15  and  18,  that  the  maximum  deviation  of  the 
electron  number  density  between  the  inviscid  core  and  the  boundary  layer 
occurs  at  the  electron- cascade  front. 


In  Fig.  B8,  variations  of  the  electron  energy  transfer  rate  thick- 
ness 6q  and  the  electron  production  rate  thickness  6ft  with  distance  x are 
shown  by  curves  A and  B,  respectively,  where 


As  shown,  these  curves  are  similar  to  each  other  and  the  largest  deviation 
between  curves  A and  B is  at  most  25$  at  Mg  = 13 .1.  Therefore,  the  approxi- 
mation of  *»  6n  in  Eq.  (3^0  is  acceptable.  6q  and  §n  are  negligible  in  the 
quasi-equilibrium  region  since  in  the  processes  of  computation  6q  and  &n  are 
written  as  (Qei  + Qinel)§Q  and  ^e^rij  0 and  - ‘ in  the  quasi -equili- 

brium region  the  ionization  rate  processes  are  in  equilibrium,  ne  (x  > xg)  -»0, 
and  consequently  the  electron  energy  transfer  rate  is  negligibly  small. 


B-3 


Wall  dissipation  terms  were  used  to  express  the  mutual  interaction 
in  analysing  the  area  averaged  quasi- one-dimensional  ionizing  shock- wave 
structure.  It  is  helpful  to  show  the  proper  expressions  of  these  terms  in 
Figs.  B9-B12. 

In  perfect-gas  boundary- layer  flows  in  a shock  tube,  the  local  skin 
friction  coefficient  Cf  is  related  to  the  Reynolds  number  by  (see  Glass  et  al. 
Ref.  Bl), 


ll-  U ?l 
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where  the  local  skin  friction  coefficient 
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and  the  shear  stress  at  the  wall 


t = 
w 


du 


For  strong  shock  waves  in  perfect  argon,  Eq.  (b8)  is  given  by 


cf  «^e  = I.91U 


(BIO) 


where  f’’(0)  = -4.062  and  Uy/Us  = 4 (see  Mirels,  Ref.  B2). 


For  a real-gas  sidewall  boundary- layer  analysis  in  ionizing  argon, 
the  shear  stress  is  given  in  boundary-layer  coordinates  as 


where 


Tw  = f"(°)  (Bll) 

W <£1 

I = j P&(x)nQ(x)Ue(x)dx 


B-U 


By  defining  the  Reynolds  number 


where 


J 

We  obtain  an  expression  similar  to  Eq.  (b8) : 


where 


As  expected,  this  equation  reduces  to  Eq.  (B8) , if  we  assume  constant  inviscid 
uniform  flow  properties  along  the  shock- tube  axis  and  C = 1 across  the  sidewall 
boundary  layer.  In  Fig.  B9,  the  local  skin  friction  coefficient  profiles  are 
shown  for  Ms  = 13  and  l6.  As  expected,  the  local  skin  friction  coefficients 
are  not  monotonically  decreasing  functions  of  x but  increase  with  x near  the 
electron- cascade  front.  Since  the  pressure  gradient  is  a maximum  there,  the 
shear  stress  increases.  Therefore,  the  value  of  Cf  v5e,  shown  in  Fig.  BIO, 
is  also  a maximum  at  the  electron-cascade  front.  This  effect  is  much 
different  from  the  results  in  the  boundary- layer  analysis  in  perfect  argon, 
and  even  greater  at  higher  shock  Mach  number  Ms  = l6,  as  drawn.  The  value  of 
cf  >jRe  at  x = 0 cm  is  about  1.58  for  both  Ms  = 13.1  and  15 -9-  Although  this  Value 
is  calculated  from  a similar  boundary-layer  analysis  for  the  frozen  uniform 
flow  conditions,  it  is  different  from  that  of  the  perfect  gas  where  cf  >jRe  = 

1.91^.  This  is  due  to  the  two- temperature  effect  of  the  argon  plasma  used  in 
the  analysis,  and  the  variation  of  transport  properties  across  the  boundary 
layer,  such  as  C = C(y)  f const,  was  taken  into  account. 

The  sidewall  heat-transfers  due  to  atom  and  electron  are  shown  in 
Figs.  Bll(a)  and  Bll(b)  for  Ms  = 13.1  and  15»9>  respectively.  The  wall  heat 
transfer  due  to  atom  temperature  gradient  -qa  is  not  a monotonically  decreasing 
function  of  x but  locally  increases  at  the  electron- cascade  front.  This  effect 
is  even  greater  at  higher  shock  Mach  number  Ms  = 16.  However,  the  heat  transfer 
due  to  the  electron  -4e  has  a maximum  at  the  electron-cascade  front. 

The  local  Nusselt  numbers  are  defined  as  follows: 


(B13) 
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where  q is  the  heat  flux  at  the  sidewall , 


q = -(x!)„  <Bll,) 

Therefore,  the  local  Nusselt  numbers  are  related  to  the  local  Reynolds  number 
defined  previously,  that  is, 


(B15) 


where  9 = Ta/Tag  for  the  atom  temperature . For  the  electron,  the  relationship 
between  Nug  and  Re  is 


Nu 


^ 0’  (0) 
a5 


(B16) 


^ ^ |1- 

. u6  ' a6  / 

where  0 = Te/Teg. 

The  ion- recombination  heat  transfer  at  the  sidewall  q^  is  defined  by 


(B17) 


The  Nusselt  number  due  to  op  and  the  local  Reynolds  number  are  related  as  follows: 


(Bl8) 


where  Z = a/ag. 

Figures  B12(a)  and  B12(b)  show  variations  of  Nu/v!Re  defined  by  Eqs. 
(B15) , (Bl6)  and  (Bl8)  for  shock  Mach  numbers  Mg  = 13-1  and  15- 9»  respectively. 
Nua/vBe  and  Nu e/*jRe  ere  maxima  at  the  electron- cascade  front.  However,  for 
the  ion-reccnbination  heat  transfer,  Nua/^e  is  a maxi  mum  at  a distance  behind 
the  electron-cascade  front.  The  reason  for  this  is  as  follows.  The  gradient 
of  the  degree  of  ionization  at  the  sidewall  Z'(0)  is  monotonically  decreasing 
with  distance  x.  However,  ctg,  at  first,  increases  from  its  frozen  value  at 
x = 0 to  the  quasi -equilibrium  peak  value  and  then  slowly  decreases  due  to  the 
radiation  energy  loss.  Therefore  OtgZ'fO)  is  a maximum  not  at  the  electron- 
cascade  front  but  at  a distance  behind  it.  This  effect  is  obvious  in  Fig. 

12(b)  for  Ms  = 15.9- 
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In  the  present  analyses,  variations  of  the  local  skin  friction 
coefficient  cf  and  the  wall  heat  transfer  -q  are  found  to  differ  greatly  from 
the  result  by  the  similar  boundary  layer  analysis  for  perfect  argon.  The 
local  Tnftvi mu  appear  at  the  electron- cascade  front  and  their  variations  are  not 
single  in  the  incubation  region  as  expected  in  the  previous  analysis.  This 
shows  the  importance  of  taking  into  account  the  mutual  interaction  terms  in 
solving  the  ionizing  shock-wave  structure  problem. 

In  the  actual  flow  field,  these  bumps  may  be  reduced  by  dissipation 
effects  but  there  is  no  experimental  evidence  that  this  indeed  happens. 
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M,  =13.1 


NON  EQUILIBRIUM 


FROZEN 


FIG.  BI(a)  VARIATIONS  OF  SIDEWALL  BOUNDARY -LAYER  DISPLACEMENT  THICKNESS  a*  AND  VELOCITY 
THICKNESS  « WITH  DISTANCE  X FOR  (a)  Ms  * 13.1,  pQ  = 5.16  TORR  AND  TQ  . 300  K, 
AND  (6)  M$  * 15.9,  p0  - 5.10  TORR  AND  T0  * 298  K.  SOLIO  LINES  ( ):  NON EQUI- 

LIBRIUM FLOW;  DASHED  LINES  (— ):  FROZEN  FLOW. 


Mj  = 15.9 


NONEQUILIBRIUM 


FROZEN 


m 


FIG.  81(b)  VARIATIONS  OF  SIDEWALL  BOUNDARY -LAVER  DISPLACEMENT  THICKNESS  6*  AND  VELOCITY 
THICKNESS  6 u WITH  DISTANCE  X FOR  (a)  M$  » 13.1,  p0  » 5.16  TORR  AND  T0  * 300  K. 

AND  (b)  Ms  • 15.9,  p0  • 5.10  TORR  ANO  T0  • 298  K.  SOLID  LINES  ( ):  NONEQUI 

LIBRIUM  FLOW;  DASHED  LINES  (— ):  FROZEN  FLOW. 


FIG.  B2  VARIATIONS  OF  SIDEWALL  BOUNDARY  LAYER  WITH  01 STANCE  X FOR  (a)  M 


13.1, 


p * 5.16  TOfiR  AM)  T_  • 300  K.  AND  (b)  M = 15.9,  P = 5.10  TORR  AND  Tq  = 298  K. 


X (cm) 

FIG.  B8  VARIATIONS  OF  SIDEWALL  BOUNDARY  LAYER  ELECTRON-ENERGY-TRANSFER-RATE  THICKNESS  6q 
ANO  ELECTRON-PRODUCTION-RATE  THICKNESS  WITH  DISTANCE  X FOR  (a)  ■ 13.1, 

p = 5. IB  TORN  ANO  T0  = 300  K,  AND  (b)  Ms  - 15.9.  pfl  ■ 5.10  TORR  ANO  T0  = 298  K. 
CURVE  A:  «<,;  CURVE  B: 


FIG.  B9  VARIATIONS  OF  SIDEWALL  0OUNOARY  LAYER  LOCAL-SKIN-FRICTION  COEFFICIENT  Cf  WITH 

DISTANCE  X FOR  (a)  Ms  » 13.1.  P0  • 5.16  TORR  ANO  T0  • 300  K.  ANO  (b)  ■=  15.9, 

p • 5. 10  TORR  ANO  T . 29B  K. 

0 o 


Nu//Re  (»/cm2) 


FIG.  B)l(b)  VARIATIONS  OF  SIDEWALL  BOUNDARY  LAYER  HEAT  TRANSFER  DUE  TO  ATOMS  qa  ANO  ELECTRONS 

qe  WITH  DISTANCE  X FOR  (a)  M$  = 13.1,  pQ  = 5.16  TORR  ANO  TQ  = 300  K,  AND  (b)  = 

15.9,  p„  = 5.10  TORR  ANO  T = 298  K. 
ro  o 


FIG.  B 1 2 (a ) VARIATIONS  OF  SIDEWALL  BOUNDARY  LAYER  LOCAL  NUSSELT  HUMBER  llu/^e  DUE  TO  qa  AND 
q AND  ION-RECOMBINATION  HEAT-TRANSFER  q^  WITH  DISTANCE  X FOR  (a)  M$  ' 13.1, 
p*  - 5.16  TORR  AND  T0  • 300  K,  AND  (b)  «s  * 15.9,  P0  * 5.10  TORR  AND  TQ  - 298  K. 


FIG.  B1 2(b)  VARIATIONS  OF  SIDEWALL  BOUNDARY  LAYER  LOCAL  NUSSELT  NUMBER  Nu/^e  DUE  TO  qa  AND 
q AND  I ON- RECOf®  I NAT  ION  HEAT-TRANSFER  qa  WITH  DISTANCE  X FOR  (a)  M$  * 13.1, 
pQ  - 5.16  TORR  ANU  TQ  * 300  K,  AND  (b)  M$  - 15.9,  pQ  * 5.10  TORR  AND  TQ  - 298  K. 
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